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We present a detailed analysis of the two-channel atom-molecule effective Hamiltonian for an 
ultracold two-component homogeneous Fermi gas interacting near a Feshbach resonance. We par- 
ticularly focus on the two-body and many-body properties of the dressed molecules in such a gas. 
An exact result for the many-body T-matrix of the two-channel theory is derived by both consid- 
ering coupled vertex equations and the functional integral methods. The field theory incorporates 
exactly the two-body physics of the Feshbach scattering by means of simple analytical formulas 
without any fitting parameters. New interesting many-body effects are discussed in the case of 
narrow resonances. We give also a description of the BEC-BCS crossover above and below Tc- The 
effects of different approximations for the selfenergy of the dressed molecules are discussed. The 
single-channel results are derived as a special limit for broad resonances. Moreover, through an 
analytic analysis of the BEC limit, the relation between the composite boson of the single-channel 
model and the dressed-molecule of the two-channel model is established. 



I. INTRODUCTION. 

Since the achievement of the first Bose-Einstein con- 
densate in a trapped dilute gas of alkah atoms in 1995 
P, 0, Q , the area of ultracold atomic gases is one of the 
most active fields of research in physics. The experi- 
mental realization of this new state of matter has been 
made possible by trapping an atomic gas cloud in a mag- 
netic trap and lowering its temperature down to about 
10-100 nK by means of a combination of laser-cooling 
and evaporative-cooling techniques. 

The extraordinary progress in the technology of trap- 
ping and cooling alkali atoms has been applied in the last 
two years also to study superfluidity in fermionic gases. 
However, in a weakly-attracting mixture of two species 
of fermionic atoms, the standard evaporative cooling 
stops to work in the degenerate regime and one can- 
not reach the superfluid transition directly. For this rea- 
son, it is more difficult to reach the fermionic analog of 
the Bose-Einstein condensation transition temperature, 
i.e., the Bardeen-Cooper-Schrieffer (BCS) temperature 
for Bose-Einstein condensation of Cooper pairs. In a BCS 
superconductor, the critical temperature of the super- 
fiuid transition is exponentially small with respect to the 
Fermi temperature, which is the temperature scale asso- 
ciated with the onset of degeneracy in the gas. Moreover, 
the BCS phase does not affect the density distribution of 
the gas in the same manner as in a trapped Bose gas. As 
a result it is not easy to detect the pair correlations of 
the Cooper pairs characterizing the superfluid properties 
of the gas. Fortunately, the use of Feshbach resonances 
in atomic Fermi gases II S S II .8, 1] , has solved il5i dfi) 
both of these problems simultaneously. 

Near a Feshbach resonance two atoms can virtu- 
ally form a long-lived bound molecular state during an 
s— wave collision. The scattering process consists of two 
incoming atoms in an energetically open channel, which 
has a different hyperfine state than the bound state in the 
closed channel. The coupling between the open and the 
closed channels is physically provided by the exchange 



interaction, i.e., the difference between the singlet 
and triplet potentials of the two alkali atoms in the 
electronic ground state. Moreover, the two channels have 
a different Zeeman shift in a magnetic field because of the 
difference in hyperfine state. Therefore, the energy differ- 
ence between the closed-channel bound molecular state 
and the two-atom continuum threshold is experimentally 
adjustable by tuning the magnetic field. As a result, the 
s— wave scattering length, and hence the magnitude and 
sign of the interatomic interactions, can be precisely con- 
trolled over a wide range. 

The possibility of tuning the interactions in a system 
is rather unusual in condensed-matter physics and has 
triggered many new experimental developments. In par- 
ticular, Feshbach resonances provide a remarkable oppor- 
tunity to study strongly-interacting fermions. By varying 
the interaction strength between the fermionic atoms, the 
crossover from a Bose-Einstein condensate of molecules 
to a BCS superfluid of loosely bou nd p air s fiol llll IT^ 
in has been explored [11 E 111 HEl IMMlIlll ■ 
These experiments have revealed evidence of fermionic 
pair condensation 0, 0] and superfluidity |23] . More- 
over, in a magnetized Fermi gas the phase separation 
between a superfluid paired component and the normal 
unpaired fermions has also been observed |H|2^. 

Feshbach resonance in atomic and nuclear physics is 
a multichannel problem. Already in 1962, Feshbach 
developed a general method to solve the multichannel 
elastic scattering equations in the presence of a bound 
state, using a formalism based on projection operators 
techniques. His approach reduces the solution of the mul- 
tichannel problem of the colliding particles to a single- 
channel problem, namely the scattering of a particle from 
an optical potential. In general, the price of this simpli- 
fication is a complicated non-local and energy-dependent 
optical potential. Neglecting the possible occupation of 
the bound state energy level and the energy dependence 
of the optical potential ultimately leads to a many-body 
theory with a local interaction described by a single ad- 
justable parameter, namely the s— wave scattering length 
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a. This is usually called the single- channel approxima- 

tion jrajg^pjiM 

IM S S IM IM EOI in the literature related to the 
study of cold atoms near Feshbach resonances. A differ- 
ent many-body approach, which allows for the population 
of the bound molecular state and thus preserves the mul- 
tichannel nature of the Feshbach problem, has also been 
considered by many authors. This is based on a two- 
channels low-energy effective Hamiltonian E3, IS 113,0, 

EllMllllSIMlEaillllEllMIIllillll where 

the effects of the resonant scattering are described by a 
bosonic molecular field coupled to the atoms in the Fermi 
sea. This bosonic field describes the bound state in the 
closed channel near the continuum of the open-channel 
scattering states. As an effect of the coupling with the 
atoms, the physical bosonic degree of freedom in the sys- 
tem is not this bare molecule but the so-called dressed 
molecule "46^ , namely the bare molecule " dressed" by the 
coupling to the continuum of scattering states. The wave- 
function of the latter is always a superposition sharing 
a component in the closed channel and one in the open 
channel. The probability of the bare closed-channel com- 
ponent is denoted by Z . Note that in this language the 
single-channel model puts Z = Q from the outset. 

In this paper we begin a systematic analysis of a two- 
channel many-body theory for a mixture of fermionic 
atoms near a Feshbach resonance based on the dressed- 
molecule picture. The paper is organized as follows. In 
Sec. ^ we introduce the physics of Feshbach resonant 
scattering in ultracold atomic gases. This preliminary 
analysis is based on a microscopic atomic physics Hamil- 
tonian with spin-dependent interactions in the presence 
of an external magnetic field. 

In Sec. Illll we introduce the low-energy effective quan- 
tum field theory above the superfluid critical temperature 
based on the atom-molecule Hamiltonian. The field the- 
ory is then solved in the so-called many-body T-matrix 
approximation. We show that this approach incorporates 
exactly the two-body scattering properties in the pres- 
ence of a Feshbach resonance. In the two-body limit, we 
find new and general analytic expressions for the bind- 
ing energy and the closed-channel component Z of the 
dressed-molecules, which reproduce experimental data 
without any fitting parameters. Next we rederive the 
many-body T-matrix approximation by means of func- 
tional integral techniques. In the limit of a very broad 
resonance the many-body T-matrix of the single-channel 
approximation is recovered. Moreover, new interesting 
many-body effects are also discussed in particular in con- 
nection with the possibility of having narrow resonances. 

In Sec. HVl and Sec. |V| the physics of the BEC-BCS 
crossover is investigated. Section llVI considers the prob- 
lem of determining the superfluid critical temperature 
approaching the transition from above, while in Sec.lVlwe 
turn our attention to the description of the crossover in 
the superfluid phase below T^. In both regimes, we give a 
systematic analysis of the crossover equations both in the 
mean-field approximation [T]| and beyond mean-field, at 



the level of the Gaussian fluctuation s IT3 . Us", "SS*! and ac- 
cording to a self-consistent approach . We focus expe- 
cially on the BEG limit of the crossover, where analytical 
calculations are possible. In the case of broad resonances, 
we show that the dressed-molecule approach |^ 0| is 
consistent with the single-channel model results. The re- 
lation between the composite boson of the single-channel 
model and the dressed-molecule of the two-channel ap- 
proach is exstablished also at the mathematical level. 

II. FESHBACH RESONANCES IN ATOMIC 
FERMI GASES. 

At low temperatures only s-wave scattering is impor- 
tant in the gas, since the centrifugal barrier prevents 
higher-order partial waves to enter the interaction region. 
Due to the Pauli principle occurring in a gas of fermionic 
particles, interaction effects can therefore only be ob- 
served when we consider a mixture of two spin states. 
For this reason we only consider such mixtures in the 
following. 

Feshbach resonances occur in atomic physics when 
studying collisions of ultracold atoms under variation of 
an external magnetic field |5^ Ig^I • In the presence of an 
external magnetic field B the spin-dependent part of the 
single-atom Hamiltonian is 

g.nt^(g)s.i + B. ^^-^-^"\ (1) 

where the energy Ohf depends on the atomic species and 
characterizes the hyperfine interaction between the elec- 
tronic spin s and the nuclear spin i. The two constants 
/ie and /ijv are respectively the electronic and nuclear 
magnetic moments of the alkali atom of interest and we 
always have that /Zg ^ ^in- 

The diagonalization of this Hamiltonian yields the hy- 
perfine states, which, at zero magnetic field, have the 
magnitude / of the total atomic spin f = s -I- i and its 
projection on the z— axis m/ as good quantum numbers. 
At high magnetic fields B 3> ahf//ie^, the hyperfine in- 
teraction can be treated as a perturbation, and, to lowest 
order, the electronic and nuclear spin projections ttt-s and 
mi are the good quantum numbers. At intermediate val- 
ues of the magnetic field, where the Feshbach resonances 
are often observed, the atomic eigenstates are a linear 
combination of the latter states, that schematically can 
be written as 

1 1 

c_i|TOs = --^]rn, = m + - > + 

1 1 

-\-ci TOc = —\mi = m > . 

2' 2 2 

This implies that only the total spin projection my is a 
good quantum number. Without much loss of generality, 
we consider a gas of ^Li atoms from now on. This species 
is most widely used in current experiments with trapped 
ultracold Fermi gases. The ^Li isotope has nuclear spin 
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Figure 1: Scattering length as a function of the magnetic field 
for two ''Li atoms in the hyperfine states |1) and |2) near the 
broad resonance at 834 G 



i = 1, which gives in increasing order of energy, the six 
hyperfine states 
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with sin0± = l/^l + (2'± + i?±)2/2, = {^in + 

2fie)B /ahf ±1/2, and = ^(Z±)2 + 2. Here the states 
\msmi) have a well-defined projection of the electron and 
nuclear spin, respectively. 

So far, we have considered only the spin degrees of free- 
dom of a single atom. The spin-dependent Hamiltonian 
of two atoms is l6l| 



Hint = (^) (si • ii 



,,.i,) + B.?^ii:^^i^, (3) 



where S = si -I- S2 and I = ii -I- i2 are the total electronic 
and nuclear spin, respectively. In view of the ongoing ex- 
periments, we consider in the following a hyperfine mix- 
ture of the states |1) and |2). Hence, the incoming chan- 
nel for two-atom scattering is the antisymmetric state 



|{12})^i=[|l)i|2)2-|2)i|l), 



(4) 



This combination can only be trapped in a far-off reso- 
nance optical trap 63] and not in a magnetic trap. It 
is, however, the most favourable mixture experimentally, 
because it cannot decay through inelastic two-body col- 
lisions. For future reference it is useful to rewrite this 



state as the linear superposition 

|{12}) =sin6l+sin6l_|l 1;1 - 1)+ (5) 
sin6l+cos6'_ ^-^|0 0;0 0) - y||0 0;2 0)^ + 
cos 6'_|_ sin 6'_ 

X (^-L|0 0;0 0) + i=|0 0;2 0)-i=|10;10) 

-I- cos 6*+ cos 0_ 1 1 - 1;1 1) 

of the basis states \SMs; lAIj), which diagonalize the 
second term in the right-hand side of Eq. Q . The quan- 
tum numbers Ms and Mj represent the total electronic 
and nuclear spin projection along the direction of the 
magnetic field, respectively. 

The central interaction for two colliding atoms depends 
on the ma gnit ude of the total electronic spin. It can be 
written as jBlj 



Vf 



Vs-Ps+Vt-Pt, 



(6) 



by using the projection operators onto states of definite 
total electronic spin S. The singlet and triplet potentials 
Vs,T are by now rather accurately known from atomic- 
physics measurements. The central interaction Vc is di- 
agonal in the basis \SMs; IMj) but induces transitions 
between different hyperfine levels. More precisely, the 
central interaction can be rewritten as 



Vc — Vdir + Vcxch, 

with the direct interaction given by 

(Vs+Vt) 

and the exchange interaction defined as 

_ (Ps-Pt)(Vs-Vt) 

* cxch — 



(7) 



(8) 



(9) 



If we describe a collision in the hyperfine basis, it is this 
latter interaction that is responsible for the coupling of 
different hyperfine states in the quantum collision of cold 
alkali atoms. Only transitions where the quantum num- 
ber Ms + Mj is conserved are allowed, since the central 
interaction cannot change the total electron or nuclear 
spin. Therefore, the state |{12}) couples only to the 
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v2 

+ sin6'+sin0_|l - 1; 1 1), 
|{36}) - -i=|0 0; 0) + -i=|0 0; 2 0) + -i=|l 0; 1 0). 

The magnetic field dependence of the s— wave scatter- 
ing length a{B) in the |{12}) channel is obtained by 
solving the system of 5 coupled-channels equations for 
these channels. The scattering length exhibits several 
resonances as an effect of the interplay between the dif- 
ferent channels in the scattering process. A very broad 
resonance occurring at 834 G is shown in Fig. ^ 

Physically, this Feshbach resonance arises when the en- 
ergy of the atoms with the incoming spin wavefunction 
I {12}), which is almost purely triplet near resonance, co- 
incides with the energy of a bound-state energy in the 
singlet potential Vs. It is convenient for this purpose to 
rewrite the interaction piece of the total Hamiltonian as 

2^e • S - ^AT • I ahf . . , • \ I -ir 

i?spin =B + (Sl • ll S2 • I2) + Vc 



Hz + -ffht 



(Si • ii -I- S2 • i2) 

+ Vc. 



-Vc 



(11) 



because \SMs;IMi), diagonalizes the Hz + Vc opera- 
tor. The hyperfine interaction _ffhf is not diagonal on 
the \SMs; IMj) states. The latter couples states with 
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Figure 2: Illustration of a Feshbach resonance. The upper 
potential curve corresponds to the closed channel interaction 
potential that contains the bound state responsible for the 
Feshbach resonance. The lower potential curve corresponds 
to the open-channel interaction potential. 



/ and iS" quantum numbers differing by one. Classi- 
cally, it gives rise to the independent single-atom pre- 
cessions of Si and U about their sum vector f^. In de- 
scribing the collision, we project the spin of the asymp- 
totic incident state |{12}) onto the new collision basis 
\SMs; I Ml) and approximate the spin state with the one 
having almost probability one. During the collision, the 
hyperfine coupling Hm can induce transitions to others, 
spin- flipped collision channels. A Feshbach resonance oc- 
curs when one of these supports a bound molecular state 
0m (r) \S'Mg; I'M'j) with energy E^, which hes near the 
continuum level of the incident |{12}) channel. The chan- 
nel of the metastable bound state is said to be closed 
because, due to energy conservation, atoms can be ob- 
served asymptotically only in the incident open channel. 
The situation is illustrated in Fig. |2 

At the high magnetic fields where the resonance occurs, 
the hyperfine interaction represents a small perturbation, 
and the total spin projections Ms and Mj are rather good 
quantum numbers. Therefore, we can approximate the 
basis with the basis \SMs]IMi) . More precisely, 

in the limit of high magnetic field we have sin 0± ~ and 
the six states in Eq. (0) can be approximated by 
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(12) 



which implies that the scattering channel |{12}) is almost 
the pure triplet state 



|{12}) 



II -1;11) 



(13) 



5 



When the magnetic field is tuned between 800 G and 
900 G, the energy of the incoming channel |{12}) lies 
near the w = 38 molecular bound state energy of the 
singlet potential Vs. As a result a very broad Feshbach 
resonance occurs. The closed channel |00; 00), supporting 
the bare molecular state, can be written by using Eq. Hl()|l 
and taking the limit of high magnetic field as 

|G0;00)^i=(|{36}) + |{25}) + |{14})). (14) 

The width of the resonance is thus determined by the 
matrix element 

(G0;00|i7M|l -1;11), 

which describes the overlap induced by the hyperfine cou- 
pling between the spin wavefunctions of the closed and 
the open channel. The total atom-molecule coupling in- 
volves of course also the overlap between the relevant 
spatial wavefunctions of the atomic potentials. 

Another resonance, a very narrow one, occurs at about 
549 G when the energy of the open channel approaches 
the energy of the same w = 38 bound state of the singlet 
potential, but now in the channel 



|00;20): 



,-l=(-LH36))-v^H25}) + -L|{14))).(l 



5) 



calculated also from Eq. (|10|) in the limit of high magnetic 
field, i.e., approximating sin0± ~ and cos9± ~ 1 at 
lowest order in Ohf //ie-B. The matrix element related to 
the width of the resonance is in this approximation 

(00;20|i7hf|l-l;ll), 

which is of the same order as the matrix element in 
Eq. H15|l . Hence, the different width of these two res- 
onances does not originate from the overlap in the spin 
parts of the wavefunctions but in the overlap of the wave- 
functions of the spatial degrees of freedom. 

In the absence of the hyperfine coupling the energy 
difference between the closed and the open channels is 
linear in the magnetic field, i.e., 



(16) 



where we have used that the difference in magnetic mo- 
ment A/i between the closed singlet channel and the 
open triplet channel is twice the Bohr's magneton or 
A/z = 2^B- The hyperfine coupling produces a shift in 
the position of the magnetic field at which the resonance 
occurs [i^ . It is this shifted value Bq that is observed ex- 
perimentally. The exact location of the resonance can be 
used to redefine the energy difference between the bound 
state in the closed channel and the threshold continuum. 
This is 



,5(5) = Am (S - Bo) , 



(17) 



and is called the detuning because it represents the vari- 
able that can be tuned experimentally by varying the 



magnetic field. At large negative detuning, the binding 
energy ein(-B) of the diatomic molecules, i.e., the energy 
you need to break up a molecule, goes linearly in the 
magnetic field B as en\{B) ~ 5{B). Near resonance, the 
hyperfine coupling leads to a quadratic magnetic-field de- 
pendence of the binding energy according to Wigner's 
formula [23 



em(S) 



ma^{B) ' 



(18) 



where the total scattering length a{B) diverges near res- 
onance as 



a{B) = Obg 1 - 



AB 

B-Bo 



(19) 



The width of the resonance AB is a phenomenological 
parameter and the scattering length a^g describes the 
nonresonant scattering in the open channel. At positive 
detuning, there no longer exists a bound state because 
the molecule decays into two free atoms due to the cou- 
pling with the atomic continuum. 

In conclusion, Feshbach resonances in atom-atom scat- 
tering involve intermediate states that are molecular 
bound states. For these molecules, the electronic and 
nuclear spins have been fiippcd from the spin of the col- 
liding atoms by virtue of the exchange interaction. Next 
we give a brief account of some recent experiments in 
ultracold Fermi gases that make use of Feshbach reso- 
nances. This account is, however, by no means complete. 
We discuss only the experimental results which are rele- 
vant in connection with the theory we want to discuss in 
the next. 

By sweeping an external magnetic field across a Fesh- 
bach resonance, it is possible to create bosonic diatomic 
molecules in an ultraco ld g as of fermionic atoms trapped 
in an optical trap [H El EUHl El III ■ The optical 
potential has the advantage to trap atoms in any spin 
state, as well as the molecules created from these atoms. 
The Fermi temperature Tp of these systems is usually of 
the order of 0.1 — IfiK. Cooling fermionic atoms down 
to temperatures T <gC Tp below quantum degeneracy is 
problematic under normal conditions, because the stan- 
dard evaporative cooling stops to work in the degenerate 
regime. However, in a mixture with two different spin 
states, the large cross section for elastic scattering near 
the Feshbach resonance can be used for efficient evap- 
orative cooling, expecially at negative scattering len gth s 
above the resonance where inelastic loss is neglig ible '71 . 
This permits experiments to reach temperatures down to 
O.OSTf. 

Trapped molecules can be created adiabatically, by 
sweeping slowly through a Feshbach resonance. If the 
temperature of the gas is comparable or lower than the 
binding energy of the molecular state, an almost pure ul- 
tracold molecular gas forms. At even lower temperatures, 
the molecules form a Bose-Einstein condensate. The pro- 
cess is reversible because the gas can be converted back 
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to atoms by reversing the sweep. These molecules are 
vibrationally highly-excited states and would usually un- 
dergo fast relaxational decay into deeply-bound states by 
molecule-molecule and atom-molecule collisions. This is 
true for Feshbach molecules made with bosonic atoms 
[22)lz3- tiiat case, the lifetime of the molecules can 
be increased only by trapping the sample in an opti- 
cal lattice 74, 75]. However, ^Li2 and '**'K2 molecules, 
consisting of fermionic atoms, show very long lifetimes 

miiEEil. 

The origin of such a long lifetime has been explained by 
Petrov et al. in 76] as follows. The size of the weakly- 
bound molecules with binding energy fv^/ma{B)'^ is of 
order a{B). This length is much larger than the charac- 
teristic radius of interaction ro . The relaxation requires 
the presence of at least three fermions at distances of 
about ro from each other. As there are only two differ- 
ent spin states, two particles are necessarily identical and 
the Pauli exclusion principle prevents this from happen- 
ing and thus suppresses the relaxation probability. 

The bound molecule has a magnetic moment which 
differs from that of the unbound atoms pairs. The differ- 
ence in magnetic moment facilities Stern-Gerlach selec- 
tion of molecules and atoms. Jochim et al. describe a 
purification scheme based on a Stern-Gerlach technique 
that efficiently removes all the atoms, while leaving all 
molecules trapped. The basic idea of this method con- 
sists in applying a magnetic field gradient B' perpendicu- 
larly to the axis of the optical dipole trap. The magnetic 
field pulls out of the trap particles for which the magnetic 
force is larger than the trapping force. The value B'^^ at 
which all the atoms are lost turns out to be smaller than 
the value B'^^^ at which molecules start to spill out of the 
the trap. This means that in principle one can remove all 
the atoms while keeping the total number of molecules 
constant. The magnetic moment of the molecules can be 
estimated from the relation /imoi = 2^at^nioi/^at7 where 
2/iat is the magnetic moment of two unbound atoms in 
the open channel. In a gas of ^Li atoms at low density 
and at high magnetic fields, the latter is almost com- 
pletely in the triplet state which implies 2/iat — 2/iB. We 
have mentioned this experiment because the magnetic 
moment of the Feshbach molecule at negative detuning 
near the resonance it will be one of the crucial quanti- 
ties in order to verify that our theory incorporates the 
two-body physics exactly. 

The long lifetime of the molecules makes it possible 
to observe a Bose-Einstein condensate of K2 jSa] 
Li2 |6E) Jfll Feshbach molecules. Zwierlein et al. cre- 
ate in Ref. [73, for example, the Bose-Einstein conden- 
sate of Li2 Feshbach molecules by operating additional 
evaporative-cooling cycles on the molecular side of the 
resonance. It is important to point out that in all these 
experiments, the condensed molecules are not detected 
directly by standard imaging techniques. In Ref. [jO] the 
molecular density is inferred from absorption images at 
zero magnetic field taken with and without dissociating 
the molecules on the positive-detuning side of the reso- 



nance. The onset of Bose-Einstein condensation is char- 
acterized by an abrupt change from a smooth distribution 
to a bimodal distribution [3, Q . Differently, in Ref. [fi^ 
Greiner et al. deduce the existence of a molecular compo- 
nent from radio-frequency spectroscopy by coupling the 
closed-channel component to a third hyperfine state ini- 
tially unoccupied. In all cases, quantitative information 
about the temperature, the total number of atoms and 
the condensate fraction, is obtained by fitting the density 
profiles using a Bose-Einstein distribution for the broad 
normal component and a Thomas-Fermi [t^ l78j] distri- 
bution for the narrow condensate component. 

The success in creating a molecular condensate de- 
pends also on a favourable ratio of coUisional rates for 
formation and decay of molecules. Atomic ^Li is optimal 
in this respect. Hence it is particularly suitable to explore 
the thermodynamics of a gas under variable interaction 
conditions. This is crucial in view of the new possibility 
of exploring the BEC-BCS crossover in ultracold atomic 
Fermi gases by tuning the magnetic field near a Feshbach 
resonance. In a dilute gas at very low temperatures the 
scattering length describes the interactions of the atoms. 
Hence, changing the scattering length by sweeping the 
magnetic field across the resonance, makes it possible to 
change the sign and the strength of interactions in the 
gas. 

When the background scattering length Obg is small 
and negative, the total scattering length is small and 
negative for large positive detuning. In this case, the 
gas undergoes a Bardeen-Cooper-Schrieffer phase transi- 
tion when lowering the temperature far below the degen- 
eracy temperature Tp. This state is characterized by a 
Bose-Einstein condensate of Cooper pairs that are largely 
delocalized in coordinate space. Their extension is much 
larger than the interparticle distance between two atoms. 

When approaching the resonance, the scattering length 
diverges and the gas enters a strong-coupling regime, de- 
fined by the condition fci?]a] > 1, where kp is the Fermi 
momentum of the gas. At resonance, the scattering 
length jumps from —00 to +00, which is the hallmark 
of a formation of a zero-energy bound state. At large 
negative detuning the superfluid gas turns into a Bose- 
Einstein condensate of tightly-bound molecules. The 
evolution between the two limits is called a crossover 
because the same U{\) symmetry is always broken and 
there is therefore no difference between the symmetries 
of the two phases. The molecular condensate obtained 
in [23, 113 represents one extreme of the BEC- 
BCS crossover. The investigation of the BCS side is 
more difficult because the fermionic superfluid transition 
hardly affects the density profile of the gas in the trap 
|26l ]. Therefore, it cannot be detected by simple imaging 
techniques. 

Several ideas have been put forward in order to 
point out a clear signature of the superfluid transition 
In m m m mil. However, it turns out that the 
Feshbach resonance itself makes it possible to study pair- 
ing effects in momentum space directly. The method, in- 
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troduced by Regal et al. at JILA in Ref. uses the 
Feshbach resonance to transfer, slowly compared to the 
two-body physics but fast compared to the many-body 
physics, the fermionic ^"^K pairs onto the bare molecu- 
lar state. This enables them to detect the condensate 
component of the gas also on the positive side of the res- 
onance. The results of the experiment of Regal et al. 
give the molecular condensate fraction after the sweep as 
a function of temperature and magnetic field. The exper- 
iments reveals a few-percent condensate fraction on the 
positive side of the resonance. The fraction depends on 
the initial temperature T of the gas and seems to have a 
threshold at about 0.5 G. This is the magnetic field at 
which kF\a{B)\ is about 1. 

The condensate component has been interpreted by 
Regal et al. as a fermionic condensate. They make use 
the absence of a two-body bound state on the positive 
side of the resonance. Hence the fermionic condensa- 
tion must be understood as a " macroscopic occupation 
of a single quantum state, in which the underlying Fermi 
statistic of the paired particles plays an essential role" . A 
similar experiment has been realized a few months later 
by Zwierlein et al. jl6| | at MIT using ^Li atoms. Surpris- 
ingly, they observe a much larger molecular condensate 
fraction than in the experiment at JILA. The origin of 
this difference still lacks an explanation. Both experi- 
ments are performed in axially elongated traps, which 
should exclude trapping effects as a main reason. 



tential /j,, and b]^ is the creation operator of a bare Fes- 
hbach molecule with momentum /ik. In addition, Vbg 
is the nonresonant or background interaction between 
the atoms in the open channel. The last two terms of 
the Hamiltonian in Eq. (|20|l describe the formation of a 
molecule from two atoms and the time-reverse process. 
The strength of this coupling is given by the bare atom- 
molecule coupling gCk). Both g and Vbg depend in prin- 
ciple also on the external magnetic field B. The energy 
<5baro = Ebarc + AfiB IS the bare detuning, where Cbare is 
the energy of a bare molecule with zero total momentum 
when the external magnetic field is absent, and A/i is the 
difference in magnetic moment between the bare molecule 
and two atoms in the states | t) and | J,). Bare quanti- 
ties must be eliminated with a renormalization procedure 
in favour of experimentally known parameters in order to 
apply this Hamiltonian to a realistic gas. In the next sec- 
tion such a procedure is introduced for the normal state 
of the atomic gas in the framework of a many-body for- 
malism based on Matsubara Green's function techniques 




A. Many-body T— matrix. 



III. ATOM-MOLECULE THEORY FOR 
FESHBACH-RESONANT INTERACTIONS. 

For an incoherent equal mixture of fermionic atoms in 
two different hyperfine states, and in the vicinity of a 
Feshbach resonance at low temperatures and densities, it 
is possible to derive [ZH 144. a low-energy effec- 
tive quantum field theory in terms of the atom-molecule 
Hamiltonian 

k.<TG{T,i} 

+ yYl ^bg(q) 4+q,T4'-q,i«k'aak,T 

k,k',q 

+ ^(y + <^barc-2M) blbk 
k 

+ ;7^J2{9* (q) &lc«t+q.i«t-q.T + • (20) 

k,q 

Here, aj^. ^ is the creation operator of an atom with mo- 
mentum h\i in the hyperfine state \a) with chemical po- 



The dynamical and thermodynamical properties of an 
atomic gas interacting near a Feshbach resonance are 
determined by the exact (four-point) interaction vertex 
function Fi^j. The latter describes the collisions of two 
atoms with spin | |) and | |), in the presence of many- 
body correlations introduced by the medium. Unfortu- 
nately, in the language of Matsubara Green's function 
techniques, the exact Fm is given by an infinite series 
of Feynman diagrams and cannot be obtained in gen- 
eral. Therefore, we have to introduce some approxima- 
tions and to restrict ourselves to the relevant subset of 
Feynman graphs. We assume that the average distance 
between the fermions is much larger than the range of 
the interactions rg, i.e., kpro <^ 1. Hence only two- 
particle collisions are expected to be important and we 
can neglect three-particle and more particle collisions. 
For a dilute system of fermions with short-range interac- 
tions we thus have to consider the ladder approximation 
HE HE 113 J which sums exactly over all two-body pro- 
cesses. In this approximation only the particle-particle 
channel is relevant and the vertex function F is reduced 
to the so-called many-body T— matrix. This satisfies a 
Bethe-Salpeter equation of the type 
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1 knT 



(21) 



h hv 



eff(ky,k,K,aOGT ( ^+k,% 



K 



T^B (k,k„K,17„) 



where Go- is the exact atomic Green's function. The ef- 
fective atom-atom interaction is given by 



Veff (kf,ki,K,f]„) 



(22) 



Vbg(k,-k.)+5*(k/)^^#^ff(k.), 



where 



Go(K,l]„) 



iMlr, 



2 



2^ - (5b 



(23) 



is the bare molecular propagator. Here ftk/ and fik^ are 
the relative momenta of the incoming and outcoming par- 
ticles of the scattering processes, ?iK is the total momen- 
tum of the two fermions, and ri„ = 2Trn/hp is the total 
bosonic Matsubara frequency of the fermion pair propa- 
gator and of the molecule. 

In general, it is very difficult to solve this equation 
directly because the interaction Vbg (k) and the atom 
molecule coupling ^(k) are included with all their details 
and are known only numerically from very complicated 
atomic physics measurements. Nevertheless, as we will 
show, in a dilute gas {kp <C ^q^) at very low tempera- 
tures, the vertex function can be calculated with- 
out explicitly solving the Bethe-Salpeter equation in all 
its complexity. This is achieved in a number of steps. 
First, we replace the complicated physical atomic poten- 
tials with a pseudopotential that leads to a more simple 
structure of the scattering equations. Then we try to 
eliminate the parameters of this pseudopotential from the 
theory in favour of an equal number of parameters that 
can be measured experimentally. The theory eventually 
works because all the relevant information to describe 
the low-energy physics of the system has been condensed 
into the fenomenological parameters. 

When atoms are colliding at ultra-low temperatures, 
the problem simplifies considerably. At these temper- 
atures they move so slowly that their kinetic energy is 
much smaller than the energy h^/mr^, related to the 
range of the interaction rg. This condition can be ex- 
pressed as 



Ath(T) » ro, 



(24) 



,1/2 



where Ath(7') = {2Tifi^ /mk-BTy " is the thermal de 
Broglie wavelength. As a result, the wave numbers of the 
relative momenta hk f and hki are much smaller than the 
characteristic wave number r^^, i.e.. 



and from scattering theory it is known that for low energy 
the T— matrix of a two-body collision becomes indepen- 
dent of the relative momenta k/,ki. Therefore, we can 
try to neglect their explicit dependence in Eq. l(?T)) by 
replacing the physical potentials with pointlike interac- 
tions Vbg^(x — x') and (?bare 5(x — x'). This implies that 
the effective interaction in Eq. H22|) can be rewritten as 



Veff (K, On) = Vbg + gbaroGo (K, On) gba 



(26) 



It turns out that under these conditions the many-body 
T— matrix is also independent of kf and ki and may be 
taken out of the integral. This gives 

T^B(K,i]„) =Veff (K,aO+ (27) 

Veff (K, n„) n (K, ao (K, o„) , 

where the kernel 



m(K,0„)=-%^^GT (^^ +k,^+c.. 



hv 



xGj, ( y ~k, — -w. 



(28) 



contains an unphysical ultraviolet divergence because the 
internal integration over k in the Bethe-Salpeter equation 
in Eq. H27|) gives a contribute for every momenta while 
in the original expression in Eq. (|21|l the natural contin- 
uum cut-off around h/r^^ of the physical potentials cuts 
off the high-momentum tail. The T— matrix renormal- 
ization procedure consists of eliminating this divergency 
by means of the the bare unknown constants (5baro, ^bg 
and (7barc from Eq. (|27() . To this purpose we introduce 
the full molecular propagator 



G (q, w„ 



Gq ^ (q, w„) - Sm (q,t^n) 



(29) 



|k/|,|k,| «ro 



(25) 



by using the molecular selfenergy defined as 

ftSni (q, ujn) = ffbarcll (q, w„) g^^ (q, w„) . (30) 

Here we have used for the momentum dependence of the 
exact atom-molecule coupling the same assumptions 
as for T'^b. For the same reason as above, there is also 
an unphysical ultraviolet divergency in the selfenergy of 
the molecular propagator that can be eliminated by a 
simple subtraction in the selfenergy. This leads to the 
renormalization of the bare detuning (5baro to the physical 
detuning 5{B) — Api{B — Bq) defined in the previous 
chapter. 
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The desired renormalized atom-molecule coupling 
is connected to the many-body T— matrix by the relation 
M 



(K, = 5bare + SbareH (K, f]„) T^^ (K, f7„), 

(31) 

where the many-body T-matrix for the background pro- 
cesses satisfies 

T^gB (K, f^„) = Vbg + ^ VbgH (K, n^) T^^ (K, n^) . 

(32) 



This means that we have dressed the bare atom-molecule 
coupling with all the two-body background interactions. 
The set of equations lf ^ - (p!^ is represented in Fig.Olby 
means of a diagrammatic notation. 



By formally solving Eq. (|2ZI) and using Eqs. ((23 EU 
and (|32|l to eliminate Go,Vbg and (jbarc, it is possible, 
by some straightforward algebra, to obtain the desired 
result 



T"B (K,aO = T^r(K,f^«) 



MB , 



(33) 



o «> 




• 
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Figure 3: The first two line define the symbol we use in our 
diagrammatic notation. The last four lines represent respec- 
tively (a) Eq. lI32J, (b) Eq. JHU, (c) Eq. (d) Eq. 



where all the bare quantities have been eliminated from 
the theory. A nice derivation of this formula by means 
of an algebra of diagrams is shown in Fig. ^ 

This result represents a generalization of the solution of 
the standard many-body T— matrix for a dilute ultracold 

§of atoms and molecules at the Feshbach resonance 
Is^ Is^ . This solution describes the scattering of a 
pair of atoms at (complex) energy iMln, while the depen- 
dence on the center-of-mass momentum K incorporates 
the Pauli-blocking effects of the medium. As we have 
anticipated, at this level the solution is still somewhat 
formal because we have only eliminated the bare cou- 
pling constants in terms of the renormalized ones. The 
next step consists in relating the renormalized quanti- 
ties, such as T^^ and g^^, to the two-body parameters. 



Figure 4: Diagrammatic calculation of the many-body T- 
matrix. The derivation starts from the definition of the many- 
body T-matrix in Eq. H27|l in terms of the bare quantities 
of the theory. The many-body T-matrix in terms of the 
renormalized couplings and propagators follows exactly from 
the algebra defined in Fig. |3] based on the exact relations in 
Eqs. ifell-ll^. 



which can be measured by experiments. 



B. Two-atom properties of the many-body theory. 

In the two-body limit [s^ |8^ [o^l there are only two 
fermions present so that in an infinite volume, the parti- 
cle density tends to zero, i.e., kp = 0, and there are no 
many-body effects anymore. This limit can be achieved 
by letting T — > and ikfln — > E while the chemical po- 
tential is fixed to ^ — 0. The dressed fermion propagator 
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Go- can be replaced by the free atomic propagator 
(k, E) ^ 



(34) 



where the energy-independent but infinite shift E^^ (0, 0) 
has renormaHzed the bare detuning (5baic to the physical 
detuning S{B) according to 



and the dependence on the total momentum K drops out 
everywhere. The many-body T— matrix in Eq. (|33|l be- 
comes the two-body T-matrix and Eq. (|33|l can be rewrit- 
ten as 



-'bare 



ftE^f (0,0) = A^B-So 



(41) 



T'^ iE)^Ti'^{E)+g'^* (E) 



G'^iE).^ 



g'^'iE), (35) 



in such a manner that the position of the resonance in the 
magnetic field is now precisely at the observed magnetic 
field value Bq. The renormalized molecular selfenergy 
can now be calculated by performing a simple integration 

M 



where the two-particle dressed molecular propagator at 
zero momentum is given by 

(E)-' = G^B (0, E)-' = Go (0, E)-' - Sf^ (0, E) . 

(36) 

The on-shell energy-dependent Tbg-matrix 



fiEf^f (k,i?)-;.s]^f (0,0) 




(42) 



/inh^ ) defines an im- 



2B f 47ryabg(g) 
J-bK [Ej = — 



1 



(37) 



is obtained by solving Eq. (|32|l in the vacuum, which 
reduces to the well-known Lippmann-Schwinger equation 

M 



where the quantity = 

portant energy scale in the problem which is related to 
the width of the Feshbach resonance. Note that as a 
result of the internal integration on k the background 
scattering length appears as an absolute value in the self- 
energy (ig. This turns out to be very important in the 
case of negative background scattering length when we 
calculate the binding energy of the molecule. The two- 
body molecular propagator takes the final form 



1 



bg 



Tl^AE). (38) hG'''{E+) ^E+-&{B)- 



|abg| 



The background scattering length abg(-B) is the quan- 
tity which is measurable experimentally and is defined 
by T2B(0) = 47rn2Qbg/m. Noting that Eqs. JSTJ and ^ 
lead to {E) = ffbaroT^^ {E) /Vbg, we conclude that 
the same kind of energy dependence can also be used for 
the atom-molecule coupling constant, i.e.. 



(43) 



Using this expression in Eq. H35|l we recover the effec- 
tive atom-atom resonant interaction in the limit of zero- 
energy scattering. This is given by 



T^B (Q) ^ 



4TTh^a{B) _ 47rft2abg(B) 2g^{B) 



(E) 



1 



1 - flbg 



(39) 



m 
m 



m S{B) 
(abg(i3)+arcs(i?)), (44) 



where g — g™' (0) is also inferred from experiments. The 
evaluation of the molecular selfenergy in the vacuum will 
complete the definition of the T^^ in terms of measurable 
quantities. 

As we have seen in the previous section, the molecular 
selfenergy in the dressed molecular propagator contains 
an unphysical ultraviolet divergency. This can be elimi- 
nated by rewriting the two-body propagator in Eq. 1361) 
as 

G^B (k, E)-^ =Go (k, E)-' + Ef,f (0, 0) 
-S^^k,£;)-E?B(o,0) 
E - ek/2 - 5{B) 



where a{B) is the total effective scattering length that 
diverges at B — Bq. 

When the background scattering length can be con- 
sidered as a constant for the magnetic field range across 
the Feshbach resonance, such as is certainly the case for 
a narrow resonance, the total effective scattering length 
can be rewritten as 



i{B) = Obg 1 



AS 

B-Bo 



(45) 



Si^k,i?)-f I]j;f (0,0) 



2B 



(40) 



with the width of the resonance AB defined through the 
relation Ai? = mg"^ / A^Th^\a\,g\A^J^■ For the broad reso- 
nance at 834G in an atomic ^Li gas, however, the back- 
ground scattering length exhibits a strong dependence 
on the external magnetic field B Wm and it is not suffi- 
ciently accurate to use Eq. H45|) i5q|. The energy of the 
molecule is given by i?k = Cm + ek/2, where the energy 
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of the zero-momentum molecular level Cm is determined 
by the pole of the retarded molecular propagator given 
in Eq. 143(1 at zero kinetic energy. This is equivalent to 
finding the zero of the equation 



Cm -S{B) ' = 0. (46) 

1 + MB)\J=%^ 



For negative detuning S{B) < the molecular propagator 
has a real and negative pole corresponding to the bound- 
state energy. The real zero of Eq. 1)46(1 can be calculated 
analytically by solving the equivalent algebraic equation. 
We find 



(7 + y/ia^ -1-72)1/3 



(7+^4^3+^2)1/3- 
21/3 



(47) 



where the coefficients are a — 3S 

3/2 



ebg + 377y/ebg and 

3 /2 

7 = — 185y/ebg — 2ei3g -f9?7ebg. The energy ebg, associated 
to the background scattering length Obg, is defined as 



Ebg 



(48) 



^bg 



Interestingly, it can be shown that the solution of Eq. 147(1 
is equivalent to the solution, in the limit when the effec- 
tive range rg can be neglected, of the equation for the 
binding energy given by Eq.(3) in Ref. j9j|. Therefore, 
in the case of a very large background scattering length, 
our theory reproduces correctly the non-trivial energy de- 
pendence due to the interplay between the resonant back- 
ground interaction and the Feshbach resonance |^ . 

In the limit of a vanishing background scattering 
length abg — > 0, the solution of Eq. ((T7|) reduces to 



SiB) 




(49) 



The bound-state energy €ni{B) of the dressed molecule 
as a function of the magnetic field is shown on the left 
sides of Figs. |5l and [7| for three different resonances. 
Near resonance the binding energy goes to zero and the 
square-root term in Eq. 146(1 , which describes the dressing 
effects of the open channel, dominates the linear term. In 
this approximation the location of the pole is at 



ma^{B) 



(50) 



according to the Wigner's formula in Eq. p8(l for the 
binding energy in resonant scattering. At large negative 
detuning, far away from resonance, the dressing effects 
of the coupling with the open channel represent only a 
small perturbation. In this limit, the energy of the Fesh- 
bach molecule approaches the binding energy of the bare 
molecule which scales linearly with the magnetic field as 



en,iB) ~ S{B). 



(51) 



In general the residue of the pole is given by 



Z{B) 



dz 



1 + 



(52) 



and is always smaller than one, because physically the 
wavefunction of the dressed molecular state near the Fes- 
hbach resonance is given by the linear superposition 



(r|xm; dressed) ~ v^Z(S)xm(r)|closed) 

, 1 e-'-/i(B) 

+ ^l-Z{B) 



y/2TTa{B) 



open) , 
(53) 



where Xm(r) denotes the wavefunction of the bare molec- 
ular state in the closed channel of the Feshbach problem. 
The dressed mo lecula r state therefore only contains with 
an amplitude a/ Z{B) the bare molecular state |xm; bare) 
of the closed channel. The remaining a/1 — Z{B) compo- 
nent is carried by the continuum of the scattering states 
in the open channel. The density of states of the dressed 
molecule at negative detuning is shown in Fig.lSJa). The 
probability Z corresponds to the wavefunction renormal- 
ization constant of the bare molecule in the closed chan- 
nel. It goes linearly in \B — BqI to zero near resonance, 
where the dressing due to the open channel is maximum. 
It goes to one monotonically at large negative detunings 
where the eigenstate of the Feshbach problem approaches 
the bare molecule of the closed channel. (See also the 
right sides of Figs.[Sl,|Bland[7I). For a broad resonance, we 
usually have that r]^ > ebg. In this case the "crossover" 
between these two regimes takes place at detunings such 
that 



Ebg. 



(54) 



When flbg is small the turning point between the 
quadratic and the linear regime is given by the condi- 
tion I Cm I — T]^- The quantity Z is also related to the 
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Magnetic Fk-ld (G) Magiictic Field (G) 



Figure 5: Binding energy and the probability Z as a function of the magnetic field near the broad resonance at 834 G for 
a mixture of ®Li atoms in the hyperfine states |1) and |2). The dashed line on the left picture represents the binding energy 
approximated according to the Wigner formula in Eq. USUI . The arrow at about 675 G on the right plot indicates the magnetic- 
field scale associated with the energy tbg ~ h^/ma^^. The curves for the binding energy and Z agree perfectly with the 
coupled-channel calculations ,6^ and the experimental data . 




Magoctif Fidd (G) Magnetic Field (G) 



Figure 6: Binding energy and the probability Z as a function of the magnetic field near the resonance at 224 G for a mixture of 
''"K atoms in the hyperfine states \ f,mf) — |9/2, —9/2) and |/, m/) = |9/2, —5/2). The dashed line on the left picture represents 
the binding energy approximated according to the Wigner formula in Eq. USUI while the stars represent the experimental data 
taken from Ref. '92*]. The arrow at about B — Bq = —1.6 G indicates the magnetic field associated with the energy fbg- The 
curves are calculated using the the same experimental parameters Obg — 174 ao, — 1.27/is and A_B = 9.7 G as in Ref. [s^. 



magnetic moment 



of the dressed molecules. Tak- 



ing the derivative on both sides respect to the magnetic 
field B in Eq. we have that 



(55) 



where we recall again that A/i is the difference in the 
magnetic moment between the atomic pair and the bare 
molecule and closed channel in absence of coupling. Us- 
ing the definition of Z, Eq. H55(l can be rewritten as 

Z(B)AAi= ^ = A^-Mmagn(S). (56) 

Therefore the change of the difference in magnetic mo- 
ments due to the dressing of the molecule near resonance 
given by the relation 



^ln^.sn{B)^A^iil-Z{B)), 



is automatically included in our theo rv l 46l| . The curve 
of Eq. (|57|l is in very good agreement j56j | with the mea- 
surement of the dressed-molecule magnetic moment in 
^Li gas taken if Ref. and discussed earlier in Sec HTl 



For positive detuning S{B) > the solution of Eq. 
has a negative imaginary part, in agreement with the fact 
that the molecule decays when its energy is above the 
two-atom continuum threshold. The molecular density 
of states at positive detuning does no longer have a delta 
function at negative energy but is given by a Lorentzian- 
like curve depicted in Fig. |HJb). The analytic expression 
which gives the maximum in the density of states em at 
positive detuning can be easily calculated in the limit of 
Obs = and reads ISJ 



1 



fj^S + 4(52 



(58) 



At large positive detuning, the peak of the Lorentzian 
distribution is located at the detuning. Near resonance, 
it is shifted from the detuning by the nonzero coupling, 
(57) and located at h^/ma^{B) exactly as the real pole at 



negative detuning. 
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Magnetic Fidd {G) 



Magnetic field (G) 



Figure 7: Binding energy and the probability Z as a function of the magnetic field near the resonance at 202.1 G for a mixture of 
""^K atoms in the hyperfine states |/, rrif) = |9/2, —9/2) and j/, nif) — |9/2, —7/2). The dashed line on the left picture represents 
the binding energy approximated according to the Wigner formula in Eq. 1501 while the crosses represent the experimental 
data taken from Ref. M. The arrow at about B — Bo = —1.2 G indicates the magnetic field associated with the energy ebg. 
The curves are calculated using the experimental parameters abg = 174 ao, — (16/9)/ifl and AB — 7.8 G as in Ref. teal. 



(s) Negative detuning 

S[E-S{B)) 

Z{B)S(E-£j^) 
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Positive detuning 

S{E-S{B)) 



Energy 



Energy i 



Figure 8: Molecular density of states in the two-body limit at negative detuning (a) and positive detuning (b). 



Many-body T— matrix revisited: functional 
integral approach. 



The path integral expression for the partition function 
of the effective atom-molecule model in Eq. ()20|l is 



The approximation developed to solve the two-body 
scattering problem, based on Eq. can be generalized 
rather straightforwardly to include the Pauli-blocking ef- 
fects on the many-body T— matrix and on the molecu- 
lar propagator. Nevertheless, for completeness we prefer 
here to rederive the many-body T— matrix approach also 
by means of functional-integral methods jOa, [93, ISg . 



(59) 



with the Euclidean action given by 



Jo J ^ 
+ dr I 0m(x, T 



V'<T(x,r) 



arc 

(B) - 2/i 



0m(x,T) 



dr 2(2m) 

+ J drjdx y dx' Vbg (x - x') Vf(x,r)V'I(x',r)i/)i(x',r)i/'T(x,r) 

/•hp p j- 

+ J dr j ds. j d-x! [g* (x - x') <;!>m ((x + x')/2, r) (x', r)7/)| (x, r) -I- c.c] 



(60) 
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where (j^m denotes the molecular bosonic field and the Grassmann fields ipcr represent the fcrmionic atoms in the two 
atomic hyperfine states. By means of the Hubbard-Stratonovich transformation 



= J d[A]d[A*] 



,+(A|'0iVT) + (V'T'>i|A) + (A|l^. i|A) 



(61) 



we can introduce the auxiliary field A(x, t), which describes the quantum fluctuations in the particle-particle channel 
according to the ladder approximation. Using this identity, we can rewrite the full action as 



hl3 



dr 2m 



V'<T(x,r) 



(62) 



+ 5 

bare 



+ J, rf-yrf-<^^(-.-)L'^ar 2(2m) 
+ dr / dx (3baro0m(x,r) - A(x,r))7/)|(x,r)?/>*(x,T) 

/■ft/3 r 

+ J dr J d:>i (pbaro <^m(x, r) - A' (x, t)) 7/)|(x, r)V'T (x, r). 



<?!'m(x,r) 



|A(x,r) 



bg 



Note that we have made the replacements Vbg(x — x') = VhgS{x — x') and g{'x — x') = .gbarc<5(x — x') as in Eq. H26|) 
above. 

Recalling the expression of the non interacting atomic Green's function 



G^,o(x,T;x',r') = jfi^ - - (5(x - x')<5(t - t') 

and using a 2 x 2 matrix (Nambu-space) notation, the dressed atomic Green's function can be written as 



Gr'(x,T;x',r') 



GfJ(x,T;x',r') 







-G^^J(x',T';x,r) 

ffbaro (/'m(x, t) - A(x, t) 

5£arcC(x,T)- A*(x,t) 

= G-i(x,T;x',r')-Sf(x,r;x',r'). 



(63) 



(64) 



X 5{x~x )S{t-t ) 



The integration over the atomic fields now involves a Gaussian integral which adds the formal result — /iTr[ln(— Gf)] 
to the effective action, that can be treated perturbatively. The effective Gaussian action is obtained by keeping terms 
up to the quadratic level in the A and fields. It reads 



S[A*,A, cj)l, =J^ drj dx A(x, 



hp 



+ I dr I dx (x,t) 

where the trace Tr(Gf,oSf)^ can be rewritten as 
Tr(Gf,oSf )^ = J^ J '^'^ J '^'^J y 

[-5barc </>m(x', t') + A(x', t')] X 

mo(x,T;x',r') [-^j^aro '/•m(x,T) + A*(x, r)] 



1 

d ft2y2 



^Tr(Gf,oSf)^ 



^dr 2(2m) 
I 



A(x, r) + 

- <5baro(B) - 2^ 



(65) 



0m (x, r), 



if we introduce the free pair propagator 

^IIq (x, t; x', r') = -G|,o (x, r; x', t') G|o (x, t; x', r') 

(67) 

in coordinate space. Using as a definition of the pairing- 
field propagator 



-;iG.Vx,T;x',r') 



5{-k-x')5{t-t') 



(68) 

no(x, r;x',r') 
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the effective Gaussian action in Eq. (|65() can be put finally in the compact form 
I 



S[A*,A,<Pl,cj,,^] ^ Jdrjcl^ j dr' j dx'{ A*(x', r') (-^G^^ Jx, r; x', r')) A(x, 



+ C(x',r') 



-ftGo ^(x,T;x',r') + -c/baro?ino(x, r;x',r')gb 



(?!)m(x,T) + 



i A*(x, T)ftno (x, r; x', r') 5baro0m(x', t') - ^(j)*^ (x, r) gbaro^no (x, r; x', r') A(x', t')|, 



r 



(69) 



which represents a generalization to the two-channel 
atom-molecule model of the Gaussian fluctuation the- 
ory developed in Ref. ^3 for a single-channel resonantly- 
interacting Fermi gas. From Eqs. (|68|) and (|32|l we see 
that the Fourier transform of the pair propagator GA^g 
can be identified with the many-body Tbg— matrix. It 
can be easily calculated by performing the sum over the 
Matsubara frequencies in XIq and using the Lippmann- 
Schwinger equation in Eq. (I38II at i? = 0. We obtain 



(70) 



^'-^Ab (-^'^")^Tbg^ (K,f2„) — 



-y 

V ^ 



K/2+K 



K/2-k 



+ 



+ 



4:TvFPahg 



The Fermi distribution factors A^k describe the effects of 
Pauli blocking in the medium. 

The renormalized pair propagator Ga, which includes 
the effects of the resonant scattering, is calculated by 
integrating out also the molecular field in Eq. 1()9|) . After 
some algebra, which involves the definitions in Eqs. H'liU 
I32|l . we recover the expected result 



Ga (K, n„) T^i^ (K, f]„) = T^^i^ (K, n„) (71) 
- g— (K, n^) — — -g- 



,MB* rx^ n \ ^ (-^J „MB j^j^ ^ ^ 



Note that in this derivation we have assumed 11 — 
XIq hich corresponds to neglecting many-body effects in 
fermionic atomic propagators. Nevertheless, this restric- 
tion would disappear if we consider the complete series 
in the expansion of the — ftTr[ln(— G)] term. At the sad- 
dle point A (x, t) — 0, there are no contributions coming 
from the /iTr[ln(— G)] term, and the solution of the pair 
propagator reduces to 



hG^' (K,r!„)= 



(72) 



-I- 



2 

3b arc 



2/i - (5baio - ck/S 



Thus the Feshbach molecule makes the pair field dy- 
namical, even before including the fluctuations around 
the saddle point, in contrast with the usual BCS theory. 



where the dynamics of the pair field is generated only via 
the fi,Tr[ln(— G)] term. As we will see in two next sec- 
tion, this can have a large effect on the bare molecular 
component above the resonance in the case of a narrow 
resonance rj^ < ep- 



D. Many-body T matrix for a very broad 
resonance. 

For a very broad resonance, the dynamics of the molec- 
ular boson close to resonance, is dominated by the self- 
energy in the molecular propagator in Eq. H29fl . The 
dynamics of the bare boson, given by the term iHu in the 
free propagator, is then effectively screened by dressing 
effects. This can be explicitly verified by carefully study- 
ing the low-energy behaviour of the exact many-body 
T-matrix in Eq. (|7TT) . 

At fixed density the width of the BEC-BCS crossover 
region, delimited by 



\kFaiB)\ ~ 



\S{B)\ 



< 1, 



(73) 



is determined by the energy scale tj'^, related to the width 
of the specific resonance under consideration. As the 
latter increases, so does the range of the magnetic fields 
which spans the crossover region. Within this range the 
detuning in general turns out, except of course very close 
to resonance, to be much larger than the other energy 
scales 



\5iB)\ » |2^|,e„ 



(74) 



To put it more precisely, at positive detuning it is 
S{B) 2/i ~ 2ei?. At negative detuning, the chemi- 
cal potential approaches half the molecular binding en- 
ergy 2\fj,\ ~ |ein| ^ l<5(-S)|. Furthermore, at low ener- 
gies fko ep and low momenta h\c[\ ~ hkp, we have 
hu! ^ hT,jn (q, ui) ~ rjH (q, ut) due to the large coupling 
rj. Under these conditions, the full molecular propagator 
in Eq. (|29|l and (|71l) can be approximated by 



ftG(q,c^)- 



/iSm (q, w) . 



(75) 



16 



(a) D^K (b) PiG-' =-S-tiY^>^<-'^CJ> 



that we obtain by integrating out the molecular field from 
the very beginning in the original action in Eq. 1160(1 . 
This shows that for a very broad resonance, the low- 
energy behaviour of the many-body T— matrix (|71|) in 
most of the BEC-BCS crossover region (I75)l . does not 
differ significantly from the many-body T— matrix of a 
single- channel interaction with scattering length a{B). A 
detailed analysis of the analytic structure of the many- 
body T— matrix H78|l has been carried out by Combescot 
in Ref. HHa^- 



X = — (o+ ° V o+ □ 



Figure 9: (a) Diagrammatic representation of the bare resonant 
interaction introduced in Eq. 1761 . (b) Approximated form of the 
molecular propagator defined in Eq. 1751 . (c) Diagrammatic deriva- 
tion of the formula in Eq. 1771 . 



Introducing also 



V,, 



(76) 



as a new definition of the bare resonant part of the po- 
tential and using the set of relations in Eqs. H26I32|I . we 

obtain without any further approximation 



(K,c^) 



^T^g^K,c.) 



g 



MB* 



(K,c^) 



1 



-<5baro - K^m (K, 
Vbg + Vrcs 
l-(Vbg+V,e«)n(K,C 



g--(K,c.) 



(77) 



This is the formal solution of a Bethe-Salpeter equa- 
tion for the many-b ody T— matrix of the bare interaction 
Vbg + Vres m il Ea. A derivation of Eq. (|77|) base d 
on the diagrammatic calculus developed in Sec. IIII Bl is 
shown in Fig. |51 Renormalizing the bare interactions to 
the two-body T— matrix as in Eq. (|7U|I this equation can 
be rewritten as 



However, even for a broad resonance, a complete ap- 
proach to investigate the physics across the Feshbach res- 
onance, cannot neglect the closed channel from the begin- 
ning. This is essentially for two reasons. First, at suffi- 
ciently negative detuning, where the systems approaches 
a gas of weakly-interacting bare molecules with binding 
energy linear in the magnetic field, the closed channel 
dominates. In that regime, the approximation introduced 
in Eq. H75|l fails because the analytic expression of the 
two-body scattering amplitude contained in Eq. H77|l has 
the wrong pole structure in the complex plane in order 
to reproduce the correct binding energy. The pole which 
denotes the binding energy of the molecular state fol- 
lows Wigner's formula e^{B) — —h'^/nia'^{B), which is 
quadratic in the detuning, at every magnetic field below 
resonance. Therefore only a two-channel model can al- 
low a full description at every magnetic field across the 
resonance. Secondly, near the resonance the gas consists 
of a mixture of atoms and dressed molecules which can in 
principle be measured separately. Dressed molecules have 
nonzero components on two different spin states config- 
urations, even above the Feshbach resonance, which can 
also be measured by experiments 20]. Thus, even when 
the dressing effect is so strong that the wavefunction con- 
tains almost completely only atomic scattering states in 
the open channel, the distinction between the two differ- 
ent components can be crucial for the interpretation of 
the experiments. 



hG^^ (K, nn)= T 



MR- 



A'!ih?a{B) 



V ^ ^ 



(K,l]„)~ (78) 

1 — AfK/2+k — AfK/2-k 1 



2/x-2ek-eK/2 2ek 



E. Many-body molecular selfenergy. 



In the limit of zero-energy scattering in the vacuum it 
reduces to the atom- atom effective interaction in Eq. (|44() 



^■Kf?a[B) 



(79) 



^{a^^(B)+aUB)) 



The physics of the dressed molecules is described by 
the molecular propagator in the presence of the medium. 
The effective Gaussian action for the molecules is ob- 
tained by integrating out the pairing field A (x, t) in 
Eq. ^ 
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'^[C'^m] ^ J dT J J dr' J dx' C(x',t') [-?iG^i(x,r;x',T' 



(80) 



1 



1 



0m (x, r). 



When we recall the definition of the pairing propagator GA^g in Eq. H68|l and Eqs. H31I32|I . we observe that the effect 
of the Gaussian integration over the field Abg consists of dressing the atom-molecule coupling constant with ladder 
diagrams. Therefore, we obtain the renormalized molecular propagator 



SI 



J dr J dxj dr' J dx' C(x',t') 



1 



nGo'(x,r;x',T') + T 3 rao(x, r; x', r') 



0n.(x,T). (81) 



By repeating the same renormalization procedure of the 
ultraviolet divergencies described above in Eqs. 1)41) |l and 
(|nj, the many-body molecular propagator can be rewrit- 
ten as 

hG^^~^{k,E) = E+2^i-S{B)-ei,/2 (82) 

where the many-body selfenergy in the ladder approxi- 
mation is given by 

fisr (q,i?)-riS2f (o,o)= 



5"^^ (2k,2ek)| 



V 



1-N. 



q/2-l-k ■ 



N. 



q/2-k 



S + 2/i - 2ek - eq/2 



(83) 

1 

2^ 



If we neglect the many-body corrections in and 
the Fermi factors in the numerator of the integrand, we 
recover the two-body result in Eq. I|42() . The integral 
of Eq. 1)83(1 cannot be performed analytically but some 
progress can be made by considering a special limit. Ne- 
glecting the many-body effects in the dressed coupling 
constant g^^, the expression reduces to 

ni:T{^,E)^hj:'^{o,o)^ (84) 



1^:15^^0,2.0 

k 



1~N. 



q/2+k 



- A^q/2-k 



£; -I- 2^ - 2ek - eq/2 



1 
2^ 



In this case, the integral can be integrated analytically 
at T = but the solution is rather cumbersome and we 
do not reproduce it here. The solution at T = in the 
limit flbg ^ 0, is given by 0| 




In 



2^ - ^ + 



^hLo + 2^j.-^ + - y^) 



2 / 



(85) 



As an important point, we must mention the fact that in 
Eq. H84|l the Galilean invariance is lost. This is because 
the Fermi sea introduces a preferred coordinate frame 
in the description of the scattering. Physical insight is 
gained by observing the T — Q solution for q = 0. In 
that case the solution is 



m:^ (o,£;)-fti]f,f (0,0) = ,? 



1 



(86) 



V 



^ , (^^V^Arctan 



- - ^^-E - 2^ Arcts 



-E-2fi 



/iSr (0, E) - fiS^f (0, 0) = 7i^-2f,-E 



(87) 



-r\\f2ep r\\J —E — 2/i Arctan 



-E-2fi 



in the limit of Obg — > when it is possible to neglect the 
energy dependence of the atom-molecule coupling con- 
stant. Only the imaginary part of the selfenergy in 
Eq. (|84|l can be calculated b y s imple analytic methods. 
In the limit Obg ^ we find [43 



V 



eiE + 2f,-^) 



X In 



coshf (^ + J(E+ + 2fi- ^)^)' 



coshf(^ 



^(£;++2^_^)^) 



The real part at nonzero temperatures can only be calcu- 
lated numerically. The effects of the nonzero temperature 
corrections are shown in Figs. 1101 and 1111 Fig. ^fa) il- 
lustrates the real part of the zero-temperature selfenergy 
in Eq. (|85|l at different momenta q. Fig. I12f b'l shows the 
effect of the corrections introduced by the background 
scattering length Obg. In Fig. IT^ a') the real part of the 
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Figure 10: Real part (left) and imaginary part (right) of the molecular selfenergy in Eq. 18411 . with q = and |(?^^ (0, 2ek) P 
\g™ (0, 0) p — g^, and at different temperatures. 





Figure 11: Real part (left) and imaginary part (right) of the molecular selfenerg y in Eq. with (0, 2ek) P 

\g^^ (0, 0) f — g^, and at different momenta and temperatures. 



selfenergy in Eq. H84|l is plotted at fixed energy ftw as a 
function of the momentum hq at different temperatures. 
Finally in Fig. ll3r b) the relative dependence between the 
two variables to and q is illustrated. 



F. Fermi edge effects for a narrow Feshbach 
resonance. 

A two-component Fermi gas with negative s— wave in- 
teractions exhibits a superconducting instability when 
lowering the temperature to T ^ Tp. The instability is 
signaled by a singularity at K = = in the two-particle 
propagator Ga which describes the Cooper channel. For 
a weak attractive background interaction Obg < 0, for ex- 
ample, the pair propagator in Eq. H7l)|l of the traditional 
BCS theory 



1 ^ r 1 - 2iVk 1 



(89) 



m 



has a pole at I l3l |9l 
Sep 



T, 



e'^ ^ exp 



2fc_F|abg| 



(90) 



where 7 = 0.5772 is Euler's constant. 

From a mathematical point of view, the logarithmic 
singularity is generated by the rather sharp Fermi surface 



in Eq. (|89|l . The same structure occurs in the selfenergy 
in Eq. (jHH- This is because many-body effects on the 
propagation of the molecule and the Cooper instability 
are described by the same ladder diagrams. As a result, 
the pairing mechanism must be enhanced by the presence 
of the molecular bound state. The effects of the resonant 
interaction on the Cooper instability, are considered in 
the next chapter, by focusing on the singularities of the 
total many-body T— matrix in Eq. (|71(l . Here we follow 
the opposite logic and concentrate on the effects of the 
logarithmic singularity, introduced by the atomic Fermi 
sea, on the molecular field. 

The bare molecular field is a true independent degree 
of freedom. Even at the lowest order in the selfenergy, 
its propagator 



hGo'{K,nn) = [ihrin- 



2m - <^ba 



(91) 



is dynamical. This is in contrast with the auxiliary pair- 
ing field Abg of BCS theory, which, at lowest order, is 
described by a static propagator 



1 



(92) 



This explains why we expect the logarithmic singular- 
ity to have nontrivial effects on the molecule thermo- 
dynamics. It turns out that it can lead to important 
consequences, expecially in the case of a narrow reso- 
nance 77^ < ep- Surprisingly, the bare molecular com- 
ponent of the gas, becomes strongly enhanced by the 
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Figure 13: (a) Real part of the molecular selfenergy in Eq. 18411 as a function of the momentum q, at different energies and 
temperatures. The various curves are (1) kT = 0, hw = ep; (2) kT = 0, hco = 0.1 • ep; (3) kT = 0, ftw = 0.01 ■ ep; (4) kT = 
0.1 • kTp, Ulu = 0.01 • ep; (5) kT = 1. ■ kTp, hcu = 0.01 • ep; (b) Real part of the molecular selfenergy in Eq. ^ as a function 
of the scaling variable hu/eq- The curves are (1) _RefiEm(t^/eq, q)/f? , (2) RefiEmii^/^cu^p) /rj , (3) -ReftEm(ti^/eq, 0.3 • ]i.p)/ri. 



presence of a sharp Fermi surface when approaching the 
resonance from positive detuning. The mechanism un- 
derlying this effect can be better understood when ex- 
amining the analogy with other well-known phenomena 
in condensed-matter physics. In particular, the resonant 
molecular level embedded in the continuum of the atoms, 
is reminescent of the Anderson model for a quantum dot, 
where, a localized electron level is located just below the 
Fermi energy of the metal leads. 

Quantum dots are small solid-state devices in which 
the numb er of electrons is a well defined integer N 
|lO0l Il0l| . In Fig. El the quantum dot is sketched as 
an electron box, separated from the leads by tunable tun- 
nel barriers and with a number of spin-degenerate energy 
states that can be single or doubly occupied by electrons 
of either spin up or down. An electron on the leads can- 
not tunnel onto the dot because it would cost an energy 
Ec — ^F, which is assumed to be much larger than the 
thermal energy k-^T of the leads. The inhibition of such 
a transition, called the Coulomb blockade, suppresses ex- 
ponentially the conduction through the dot at low tem- 
perature. This suppression occurs because the process of 
electron transport through the dot involves a real tran- 
sition to the state in which the charge of the dot differs 
by one unit from the thermodynamically most probable 
value. However, the Heisenberg uncertainty principle al- 
lows higher-order processes for short durations, in which 



(a) 



SE\ 



(b) 




Figure 14: Examples of co-tunneling processes: (a) inelas- 
tic co-tunneling transferring an electron between the leads 
leaves behind an electron-hole pair in the dot; (b) elastic co- 
tunneling; (c) elastic co-tunneling with a flip of the spin. This 
figure is taken from Ref. |101| 



virtual states participate in the tunneling process. The 
leading contributions to this activationless transport are 
provided by the inelastic and elastic co-tunneling pro- 
cesses described in Fig. E] 

It turns out that in a quantum dot the amplitude of the 
elastic co-tunneling process with a singly-occupied level 
of Fig El (c), diverges logarithmically when the energy 
hw c^i k^T oi an incoming electron approaches zero. The 
singularity in the transition amplitude gives a dramatic 
enhancement of the conductance through the dot. In the 
elastic co-tunneling process (c) the electron on the dot is 
quickly replaced by another electron, when the electron 
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on the dot tunnels to one of the leads. The characteristic 
time scale for such a co-tunneling is about h/Ec- Events 
of type (c) also effectively flip the spin on the dot. Suc- 
cessive spin-flip processes screen the local spin on the dot 
until a spin singlet is formed by the electrons in the leads 
and on the dot. Macroscopically, the system enters in a 
many-body correlated quantum state. The formation of 
this entangled state, which is a pure many- b ody effect , 
represents the Kondo effect in quantum dots |lO0l Il02| . 
If we interprete the tunneling as an effective magnetic- 
exchange coupling, the physics of a quantum dot between 
two leads becomes analogous to the original Kondo effect 
|l03j| for magnetic impurities coupled to conduction elec- 
trons in a metal host. As shown in Fig.^J^a), the Kondo 
effect is signaled by a narrow resonance at the Fermi en- 
ergy in the density of states of the dot. 

Returning back to our model, we expect something 
similar to happen when the molecular level is located at 
somewhat more than twice the Fermi energy above the 
threshold of the atomic continuum. Differently than in 
quantum dots, the molecular state has to lie above twice 
the Fermi energy of the atoms, because otherwise the 
ground state of the gas would contain a Bose-Einstein 
condensate of molecules. The atom-molecule coupling 
produces virtual transitions between the molecules and 
two atoms above and two holes below the Fermi sea, as it 
is shown in Fig. 1161 In this case the virtual co-tunneling 
of pairs of atoms to the molecular level is indeed remi- 
niscent of the inelastic co-tunneling of Fig.E|(c). How- 
ever, in our case, no Coulomb blockade is required. The 
Fermi sea removes the symmetry between scattering to 
and from the virtual molecular level, because atoms can 
tunnel only in pairs. This leads to the logarithmic singu- 
larity in the selfenergy in Eq. (|84|) as mentioned above. 
The onset of a singularity in the selfenergy signales the 
formation of a new many-body collective state. Taking 
the lowest order precesses shown in Fig. 1161 this many- 
body resonant state can be written as a linear superpo- 
sition 



|k|>*:F 



/3k&o4«k,Ta-k,il*F) 



|k|<feF 



, (93) 



where |$f) indicates the filled Fermi sea and the molec- 
ular field is treated in the zero-mode approximation that 
considers only the lowest molecular state. This is justi- 
fied at very low temperatures. In the formation of the 
many-body coherent state the system gains an energy 

^Ek ^ {^k\H\-^k) - Eg. (94) 

where Eq = 2X)|k|<fcF(^k — m) is the ground state of the 
Fermi mixture in absence of Feshbac h res onace. This 
shift can be calculated variationally |l04i ^05] as fol- 
lows. First we evaluate iJl^x) considering the Hamil- 
tonian given in Eq. I|102|) . Then we project this state 



on the three different components l^'K.j) of the ground- 
state ansatz of Eq. (|93|l . As a result Eq. 1)94(1 is turned 
in three equations for the variables A£', at, /?k- The 
variable Z cancels out because it is the normalization 
constant. Eliminating ak and /?k we obtain at T = a 
self-consistent equation for Ai?. This is 
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2/i - 5{B) 



(95) 



i]-^/I^Ek 



y/AEK+2lI+^' 



where we have used that ji — ep = h^kp/2m. At large 
detunings 5 ^ 2)1 this equation has only one solution 
at /S.Ek — 5{B) as it is expected when resonance ef- 
fects can be neglected. However, when 5 > 2/i, we have 
three solutions. The lowest energy solution defines the 
ground state of the systems and occurs when I^Ek — 0. 
The logarithm in the last term diverges when AEk 0. 
Therefore, we have that 2fj. > AEk, \ In AiJ^I > AEk 
and we can approximate Eq. H95|) as 



2n-S{B) ~?7 



2y/2jl 



2 + In 



\AE 



8n 



(96) 



This leads to a non-perturbative expression for the energy 
shift of the ground state given by 



\AE. 



8n 



(97) 



The change in the ground state energy defines a temper- 
ature, i.e., the analogous of the Kondo temperature. 

The singular behaviour of the selfenergy is associated 
with the occurence of a resonance in the spectral density 
p(q, ui) of the molecular bosons given by 



>(q,c.) = -Im[G(+)(q,c.)]M. 



(98) 



This feature is clearly illustrated in Fig. llSf bl. 
Fig. El shows the function /9ni(q, w) for different detun- 
ings and momenta q both at nonzero and zero temper- 
ature. At large positive detuning, the spectral function 
shows just a single broad peak centered around the de- 
tuning. This is the expected situation for a single molecu- 
lar state with a finite lifetime. As the detuning gets closer 
to twice the Fermi energy, the spectral density shows two 
other sharp peaks slightly above and below the zero fre- 
quency. When approaching the resonance, the spectral 
density of the molecular field at low momenta and tem- 
peratures does no longer satisfy the unitarity sum rule 



d{hhj) p^{q,uj) 



1, 



(99) 



despite the fact that the selfenergy in Eq. H84() is analytic 
in the upper-half complex plane. This happens already 
above the BCS critical temperature and is due to the 
appearance of a complex pole in the upper half plane in 
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Figure 15: (a) Kondo resonance in quantum dots; (b) Many-body resonance in the molecular density of states near the Feshbach 
resonance 




Figure 16: Virtual tunneling of pairs of atoms to the molec- 
ular level in the Fermi gas near a Feshbach resonance. 



peak. 

The additional spectral weight at low frequencies, in- 
duced by the above mentioned many-body effects, leads 
to an increase of the bare molecular component in the 
gas, calculated by means of 

nB(T) = /"J d(M / (0iP4q,-)iV (^) . 

(100) 

with respect to that estimated from two-body physics. 
The enhancement of the bare molecular components for 
realistic resonance parameters is discussed in detail in 
Ref. 113. 



the propagator. This feature is an artifact of the weak- 
coupling approximation, which neglects self-consistency 
effects on the Fermi propagators in the calculation of 
the molecular selfenergy when the gas enters the strong- 
coupling regim e kpa{B) < —1. For fermions, Luttinger's 
theorem jlOTj shows that the causality of the selfenergy 
is a sufficient condition to have also a causal propaga- 
tor. We are not aware of an analogous theorem in the 
case of b osons and our findings indeed present a counter 
example Il08|. The double-peak feature arises because 
the spectral density must be negative at negative fre- 
quencies for a bosonic field. By plotting the quantity 
/3m(w)/w |l09j |. as shown in Fig. ITHT a). a single peak 
at zero frequency is obtained. We call this a molecu- 
lar Kondo resonance. The effects of nonzero momentum 
and of corrections due to a nonzero background scatter- 
ing lenght flbg are shown in Fig. llSt b). Both tend to 
suppress the resonance. Nonzero momentum corrections 
spoil the coherence of the many-body virtual cotunneling 
events similarly as BCS pairing is suppressed at nonzero 
total momentum. A large background scattering length 
means that there is a bound state in the open channel po- 
tential ■ This induces some density of states at large 
energy. Because the total density of states is normalized 
to 1, this results in an effective reduction of the resonance 



G. Exact mapping to the anisotropic Kondo 
Hamiltonian. 

The analogy with the Kondo effect in a quantum dot, 
introduced in the previous section, requires a more for- 
mal analysis. It is very desirable to search for an ex- 
act mapping from the atom-molecule Hamiltonian to the 
Hamiltonian of the Kondo problem. 

The physics of the Kondo effect is described by a model 
Hamiltonian based on the assumption that the magnetic 
moment of a local impurity is coupled via an antifer- 
romagnetic exchange interaction J with the conduction 
electrons. This is known as the Kondo Hamiltonian and 
reads 

H^Y. (£k-M)<.Ck,. [j+5+c+^Ck',T (101) 

k,cr k,k' 

+J_S'"c+^Ck'a + JzS'^(cJ|Ck',T-c+j^Ck/^j^) , 

where ek is the energy of the conduction electrons and 
Sz and >S'*(= Sx ± iSy) are the spin operators for the 
impurity with spin 1/2. The operators Ck.a and c^^ 
correspond to the creation and annihilation operators of 
conduction electrons with momenta k and one of the pos- 
sible spin states a = | t), or ja =|) which scatters of the 
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Figure 17: Molecular spectral density (left) calculated from the selfenergy in Eq. H87|l at different detunings. Molecular spectral 
density (right) calculated from the selfenergy in Eq. H84^ at 5 = 5e_F for different temperatures; (1) T = Tf, (2) T = O.STf, 
(3) T = O.ITf, (4) T = 0.06Tf, (5) Tbcs ^ OmkTp. 
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Figure 18: (a) The spectral density is plotted together with p{uj)/u. The double structure of the molecular resonance merges 
then into a single Kondo-like resonance, (b) Effects on the molecular spectral density, calculated from the selfenergy in 
Eq. (1871 . of the kinetic energy of the molecules: (1) g = 0, (2) g ~ 0.3 ■ kp and (3) q — kp- The curve (4) is calculated from 



the selfenergy l|86|l using q = 0, 
background scattering length. 



dbg ~ 10 ao, ep ~ 3.6fiK. It shows the corrections given on the density of states by a large 



impurity. The coupling constants J± and Jz describe 
spin-flip scattering. 

The Hamiltonian of our model in Eq. (|20|) appears 
rather different than the Hamiltonian in Eq. (|101|l . but it 
becomes more similar vi^hen a restricted many-body Fock 
space is considered ^110]. For our present purpose, we can 
in first in stanc e consider only the lowest-energy molecu- 
lar state We can also neglect the efl'ects of the 
background interaction. Within these approximations, 
the atom-molecule Hamiltonian in Eq. (|20ll becomes 

H= (ek-M)4,aak,^ + ('^(S)-2/i)6j6o+ 

k,<Te{T,i} 

+ E(3* ^oak,ia-k,T + h.c.) . (102) 

k 

The bosonic number operator bj&o s-cts on the occupied 
(unoccupied) lowest energy molecular state |1) (|0)) ac- 
cording to b^boll) — 1|1) and b^bolO) = 0. Hence, in 
this reduced Hilbert space, we can formally identify the 
molecular operators b\ and bg with spin 1/2 raising and 



lowering operators 5+ = 5^ + iSv and S- ^ S"" - iSy , 
which satisfy analogous relations 

S+S-\ + 1/2) = Q + I + 1/2) = 1| + 1/2) (103) 

S'+S'-| - 1/2) = Q + 5-"^ I - 1/2) = (104) 

in the two-dimensional Hilbert space {| + 1/2),| — 1/2)}. 

More progress can be made by rewriting the atom- 
molecular coupling in terms of the symmetrized atom- 
pair operators dk ~ a-k,iak,| + ak,j.a-k,T/V2 and = 

'^k t'^-k i+'^-k J"^i which satisfy the commutation 
relation 

[<, rfk] = \ (+iVk.T + ^-k,i + iV_k,T + N^.X - 2) . 

(105) 

In the restricted two-dimensional Hilbert space, defined 
by the condition that the eigenvalues of A^k.T + ^-k.i 
and -/V_k,| -I- -^Vk,| are either 2 or 0, the commutator in 
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Figure 19: Scaling trajectories for the anisotropic Kondo model 
in Eq. tToTt [TT^ . 



Eq. (|105|) can be put into the matrix form [d^ , 
a^, by means of the z-component Pauh matrix 



1 

-1 



(106) 



This suggests the correspondence ^ = |(o'^+*cr|^) 
and ^ cjy^ = ^{(J^ ~ *^k)- Finally, by rewriting the 
spin ladder opera tors , by means of the Abrikosov rep- 
resentation |ll4j . in terms of anticommuting fcrmionic 
operators = cj^|Ck,x and a^^ = cj^j^Ck,| the atom- 
molecule Hamiltonian in Eq. (|2()|l can finally be trans- 
formed into 

= 51 - <„Ck,a + (<5 - 2/i) (107) 
^[S'+c+^Ck,i + S'-c+^Ck,T]. 



When S — 2fi, this Hamiltonian is equivalent of the 
anisotropic Kondo Hamiltonian in Eq. p01|) if we state 
the correspondence g/V^ — J+ = J- and Jz = 0. At 
finite detuning but larger than 2/i, the term (S — 2/i) 5*2 
represents an external magnetic field applied on the local- 
ized spin, which removes its degeneracy and ultimately 
destroys the Kondo effect. Moreover, it is important to 
realize that, the lack of the nonflipping spin- interaction 
term, signaled by the condition Jz = 0, does not rule out 
the possibility that such a term is generated by higher- 
order processes. As illustrated in the example of Fig. 1191 
the renormalization group flow in the solution of the 
Kondo problem, can generate n onze ro values of Jz start- 
ing from a point with Jz = 'll3]. A renormalization 
group analysis of the coupling constants flow reveals that 
this occurs also in Eq. H107() . 



H. Fermi liquid description of the normal state in 
the weak coupling regime. 

In absence of the resonance, the ultracold dilute 
fermionic mixture weakly-interacting in the background 



scattering channel above some certain critical tempera- 
ture Tc, is a Fermi liquid |ll5l |. In this section, we con- 
sider the weak coupling regime l/kFa{B) < —1 in pres- 
ence of the Feshbach resonance from the point of view of 
the Fermi liquid theory ,49j . The thermodynamics of a 
normal-fluid Fermi liquid is determined by the fermionic 
Green's function 



ihWn - (fk - - ftS^ (k, LUn) 



(108) 



In a self-consistent approach, the fermionic selfenergy 
hT,l{'k,uJn) is calculated in the ladder approximation 
from 



^I^^(k,C.„) = -^/d3k'I- 



E 



(109) 



xTMB(k + k',u;„-l-0 G_,(k',0, 

where the many-body T-matrix of Eq. (|33|l is deter- 
mined by Eqs. 1(221), and ^7^. From the selfenergy 
?il]^(k, a;„) we can calculate the atomic spectrum given 
by the p oles of the fcrmionic propagator in Eq. H108|) 
[Sa. IllSj according to the equation 



Wk 



2m 



ftE^(k,cjk) = ek + ft5]-''(k,t^k), (HO) 



where we have suppressed the spin indices because the 
two-species have the same mass and chemical potential. 
This equation can be simplified in the weak coupling limit 

as 



^k ^ Ek + ftS-^(k, ek), 



(111) 



because we expect a quasiparticle spectrum for the single 
particle excitations of the type Wk = ek[l + 0{kFa{B))\ 
and the selfenergy ?iS-^(k, Wk) to be of the order of 
0[kFa{B)). The real part Re[;il]^(k, ek)] gives a shift 
in the single particle energies of momentum k, while the 
immaginary part Im[wk] = ft7k leads to a finite lifetime 
proportional to l/7k- Since the lifetime becomes infinite 
at at kp^ the chemical potential is defined by 



hJ:^{kF,ekF)- 



(112) 



Close to the Fermi surface the energy spectrum can be 
expanded into a Taylor series 



Re[t^k] = ffcF + -{k - kp) + 0{k - kpf, (113) 

m* 

which defines the effective mass 



m = n kp 



\ dk \k=ki 



(114) 



The effective mass determines the properties of the sys- 
tem in the zero-temperature limit. 

In presence of the resonance it is important to dis- 
tinguish |49| between the weak-coupling regime near a 
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broad and a narrow resonance. For a broad resonance 
the condition ^ ep implies that in the weak-couphng 
regime l/kpa{B) < —1 we have S{B) ^ ep- Therefore, 
we can neglect the retardation effects on the atom-atom 
interaction induced by the dynamics of the bare molecu- 



lar boson. According to the discussion of Sec. lIIIDl near 
a broad resonance, the many-body T-matrix of Eq. 
can be approximated by the single-channel expression of 
Eq. H78|) . Then the atomic selfenergy becomes 



{2n) 



-y 



m 



a{B) 



47rft2 1 

TO V 



E 



1 - N, 



(k+k')/2-l-p 



N, 



(k+k')/2-p 



iujn + 2// - 2ep - e(k+k')/2 



1 

2Z 



G.(k',c.;) 

(115) 



In the weak-coupling limit the many-body T-matrix in 
Eq. (|115|l can be expanded in powers of a{B) and we 
obtain 



n ATih^a{B) 

2 TO, 



(116) 



where n is the total atomic density of the gas and 
/iS^(k, a;„) represents a correction of order {kpa{B))^ . 
Calculating the selfenergy hT,^ {k, ojn) in the limit of zero 
temperature one finds the expression of Wk, 7k and m* 
up to order {kpa{B)Y of a dilute Fermi liquid theory 



TO 

m 



{kpa{B)) 



2m ' Svr 
-^{U-2ln2)(kpa{B)f 



= 1 



7k = 



hkl 2 

2m TT 



{7\n-l){kFa{B)f 

kp — k^ 
kp 



(117) 



for \k - kp\/kp < 1. 

In the case of a narrow resonance the bare dynamics 
of the molecular propagator cannot be neglected in the 
weak-coupling limit \/kpa{B) < —1. Therefore, we have 
to retain the full expression of the molecular propagator 
in the many-body T-matrix of Eq. (133(1 . Using the formal 
definition of the effective mass in Eq. pi4|l and calculat- 
ing the integral in Eq. (|109|l . one recovers the additional 
contribution to the effective mass of Eq. 1117(1 



5m* 



9 



2 S{By 



(118) 



calculated by Bruun and Pethick in Ref. H^. How- 
ever, we have seen that, in the weak-coupling regime 
l/kpa{B) < —1 on the positive side of the resonance, 
non-trivial correlations arise in the molecular degree of 
freedom. These correlations induce a resonance in the 
molecular spectral function. As a result the finite-lifetime 
molecule on the positive side of the resonance cannot 



be pictured as a well-defined quasiparticle as one would 
naively think in the case of a narrow resonance. There- 
fore, it is reasonable to ask if it is not possible that the 
Fcrmi-liquid description could partially break down in 
the case of a narrow resonance in combination with the 
onset of the molecular resonance at the Fermi surface. 
This requires a more detailed analysis of the fermionic 
Green's function in Eq. 1(11)8(1 and will be addressed else- 
where. 



IV. THE EQUATION FOR THE CRITICAL 
TEMPERATURE Tc IN THE BEC-BCS 
CROSSOVER IN RESONANT ATOMIC FERMI 
GASES. 

The ability to control experimentally the strength of 
the interactions between ultracold atoms, offers the excit- 
ing possibility to study in detail the crossover between the 
Bose-Einstein condensation (EEC) of diatomic molecules 
and the Bardeen-Cooper-Schrieffer (BCS) transition of 
atomic pairs [Ullllllllllillillllig. 



A. Thouless criterion. 

According to the Thouless criterion |ll6l |. the super- 
fluid phase transition in a Fermi system with attrac- 
tive interactions, occurs at the temperature for which 
the many-body T-matrix in Eq. I(33() develops a pole at 
K = cj = 0. At large positive detuning, i.e., in the 
weak-coupling BCS limit, the Thouless criterion yields 
the equation for the determination of the BCS critical 
temperature of a Fermi gas with attractive interactions. 
However, when approaching the resonance, the attraction 
increases and the Fermi surface of the gas gets strongly 
renormalized. Hence, the approximation /i ~ e^, valid 
in the weak-coupling BCS limit, is no longer valid. Par- 
ticle number conservation gives the desired condition on 
the chemical potential of the gas under equilibrium con- 
ditions 0, 0, 01 . This leads to a set of two coupled 
equations for Tc and fi. In the opposite BEC limit the 
role of the two equations is inverted. The equation of 
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state determines the critical temperature for the Bose- 
Einstein condensation of dressed molecules. The pole of 
the many-body T-matrix becomes two-body in nature. 
This fixes the chemical potential to half the binding en- 
ergy of a dressed molecule. 

The equation of state can be determined in differ- 
ent approximations. The simplest mean- field theory at 
nonzero temperature, neglects finite-lifetime effects of the 
molecules and does not reproduce the correct binding en- 
ergy of the dressed molecules. The latter problem can 
be repaired by replacing in the equation of state the 
bare energy of the molecular level with its dressed value 
calculated at Sec. IIII Bl in the two-body scattering ap- 
proximation. This leads to the modified mean-field de- 
scription 'isj discusssed in section HV Bl However, this 
mean-field approach still neglects finite-lifetime effects of 
the molecules. As a result, it exhibits a phase transi- 
tion between the BCS phase and the BEC phase. Since 
we expect, from symmetry principles, a smooth curve 
for Tc and fi, this must be considered to be an arte- 
fact of the mean-field theory. The inclusion of Gaussian 
fluctuations on top of the mean-field solution, according 
to the Nozieres-Schmitt-Rink approach [13 . 11^ . is suf- 
ficient to obtain a smooth crossover in the two-channel 
model 4^ l47l l5fl| | . This approximation is described in 
Sec. IIV CI A comparison of the critical temperature as a 
function of the magnetic field for a broad resonance re- 
veals that our modified mean-field approach agrees quan- 
titatively fairly well with the Nozieres-Schmitt-Rink re- 
sult. The Gaussian approximation exhibits a maximum 
in the curve for Tc |45l. ISOl IllTj . It has been shown that, 
in a self-consistent approach, this maximum in the curve 
for Tc disappears both in a single-channel 0| and in a 
two channels model [s^. In Sec. IIVDI we discuss the 
self-consistent equations in connection with our dressed- 
molecule picture through a systematic analysis of the 
BEC limit where analytical calculation are possible. 



B. Modified mean-field equations. 

In order to study the pole structure of the many-body 
T-matrix in Eq. H33I) . it is crucial to consider explicitly 
the many-body corrections of the atom-molccule coupling 
constant in the calculation of the selfenergy of the dressed 
molecular propagator in Eq. H83|l . However, a more trans- 
parent derivation of the equation for the critical tem- 
perature can be obtained by rewriting the many-body 
T-matrix in Eq. H33|l in terms of the simpler two-body 
T-matrix in Eq. (|35|l . We start by considering the defi- 
nition of the many-body T-matrix in Eq. I|27|) 

T^^ (K, E) = Voff + VoffH (K, E) T^^ (K, E) , 



with Vcff given in Eq. H26() . Then we replace the kernel 
in Eq. ^ 



hv 

with the result of Eq. ^ 
kBT 
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E 



(119) 



nuo iK,E) 



hV 



K 



E 



X Go,i — - k, — 



2 ' 2 

S + 2/1 - 2ek - eK/2 



(120) 



K/2-k 



as we have done in the calculation of the molecular self- 
energy in Eq. (|83|l . Subtracting from the many-body 
T-matrix equation in Eq. (|27l) the two-body limit 



T^^ (i?) = Veff + Ve^ E 



1 



E-2e^ 



T^^ (E) (121) 



at E = 2/Lt, the many-body T-matrix can be put into the 
form 



T*^^ {K,E) 



(122) 



T^B {2^l+) 



l-T2B(2/,+ ) no(K,i?)-i.Ei 



2u+-2ek 



where the two-body T-matrix T^^ {E) is known from 
Eqs. (|35I43|) . At K = w = this becomes 



T^B (2/i+) 



1 + T2B (2/i+)S(0,0+)^ 



where 



K/2+k 



K/2-k 



£; + 2/1 - 2ek ~ eK/2 



(123) 



(124) 



According to the Thouless criterion the denominator of 
Eq. H123|) vanishes at the critical temperature. To illus- 
trate more clearly the physics involved we consider first 
a system when the background interaction can be ne- 
glected, i.e., the limit Obg 0. This is the case, for 
example, in '*'^K mixtures where the background scatter- 
ing length is small and does not contribute to the pairing 
mechanism because is positive. In that case the condition 
to have a pole in Eq. 11231) becomes 



SiB) - 2/i + 77 - (0, 0+) = 0. 



(125) 



The critical temperature Tc and the chemical potential 
fi are determined at every magnetic field in the mean- 
field approximation by solving self-consistently Eq. H125|) 
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together with the equation of state. At the mean-field 
level, neglecting the finite lifetime of the molecules, the 
latter is given by 



"tot 



1 dfl 



(126) 



y [2A^(ek - /i) + 27VB(ek/2 + 5{B) - 2/i)] 



where the free-particle atomic and molecular propaga- 
tors are given by Go/h = {iHuj — 6^^/2 + 2^ — d{B))^^ and 
Gafi/h = {ihijj — Ck + Ai)~^i and Nb is the Bose distribu- 
tion function. However, to neglect completely selfenergy 
effects in the molecular propagator is not consistent. This 
can be repaired in a first approximation, by replacing the 
detuning 5 in the free molecular propagator occurring in 
Eq. H126|) with the two-body curve of Eq. (|S5|l 




(127) 



that we have discussed in Sec. 2.3. This replacement is 
necessary because the energy of the dressed molecule in 
the equation of state must be the same as the binding 
energy resulting from the chemical potential in the BEC 
limit from the gap equation. As we shall see, this is cru- 
cial in the case of a crossover near a broad resonance 
where the binding energy is quadratic in the detuning 
in a large part of the crossover region. For a very nar- 
row resonance, when 77^ <C e_F , the binding energy of the 
dressed molecule can be approximated as em ~ 5{B) ex- 
cept that on a vanishing small range of magnetic field 
close to resonance and the equation of state in Eq. H126|) 
turns out to be a good approximation. However, the res- 
onances used in current experiments on the BEC-BCS 
crossover are broad in nature. Therefore, the mean-field 
equation of state has to be modified to 

n = ^ [2iV(ek + 2NB{e^/2 + - 2^^)] . (128) 

k 

Note that the two-body quantity em does not account 
for a possible many-body shift of the molecular bound- 
state. However, for a broad resonance, this seems to be a 
smaller correction when compared to the two-body shift. 

At positive detuning, in the BCS limit, the chemical 
potential is positive and the many-body term of Eq. H125() 



g^E (0, 0+) = g^Re [S (O, 0+)] + 77^^. (129) 
Therefore we obtain the gap equation 

S{B) -2 fi^ g^Re[E {0,0+)] . (130) 



Note that our derivation is based on the many-body T- 
matrix calculated in the ladder approximation. There- 
fore, it neglects the Gorkov's corrections to the tran- 
sition temperature introduced by the effects of the den- 
sity fluctuations of the medium on the effective two-body 
interactions. 

It turns out that for a narrow resonance rf ~ ep, we 
can be in the weak-coupling regime kp\a{B)\ < 1 at mag- 
netic fields such that (5 ~ 2/i ~ e^?. In this case the BCS 
critical temperature is given by 



8eF -v-2 2.:^ 
e ' e " "I 



exp 



2fcF|a(-B)| 



(131) 



When the resonance is broad and 77^ ^ 2eF, the detun- 
ing is always (5 ^ 2/x ~ ep in the weak coupling range 

kp\a{By\ < 1, and the term c 1 can be neglected. In 
this limit the BCS critical temperature can be approxi- 
mated by 



e ' exp 



2kF\a{B)\ 



(132) 



This is the same result that one would obtain by looking 
at the pole of the the single- channel many-hody T-matrix 
given by Eq. (|77|l . Therefore, we conclude that the posi- 
tive shift of the weak-coupling BCS critical temperature 

by the factor e " >' , characteristic of the resonant super- 
fiuidity, can be observed only in narrow resonances. 

At large negative detuning, the pole of the many-body 
T-matrix in Eq. H27|l determines the chemical potential 
2fi of the dressed molecule, i.e., the energy we need to 
make a dressed molecule from two atoms at zero mo- 
menta. In the BEC limit this is very large and negative 
and corresponds to the binding energy of the dressed 
molecule. Indeed the pole of the many-body T-matrix 
becomes two-body in nature because, for negative /i, the 
many-body factor is 



g^E (0,0+) = .g2Re[S (0,0+)] 



(133) 



and vanishes exponentially as \2fj,\/k-BT — > 00. The con- 
dition (|125|) becomes 



(5(B) - 2^ 77^^ - 0, 



(134) 



which is the same algebraic equation as in Eq. 146() , again 
in the limit of Obg 0, satisfied by the binding energy 
in Eq. 



(135) 



Notice that the curve em of Eq. p27|l . that we have re- 
placed in the equation of state coincides with the binding 
energy of Eq. H135II at every negative detuning. 

Once more, it is instructive to distinguish between the 
two limits of narrow 77^ <C €p and broad ep <C 77^ reso- 
nances. On the BEC side of a narrow resonance, when 
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the chemical potential approaches half the binding en- 
ergy we have n ~ em/2 ~ (5/2, and the dressed molecu- 
lar wavefunction has only a small amplitude in the open 
channel, that is Z ~ 1. This can be checked easily if we 
assume that the chemical potential approaches half the 
molecular binding energy /i ~ em/2 approximatively at 
some value -Bbec when fci?|a(SBEc)| — 1- From the def- 
inition of the resonant scattering length in Eq. (|44|l . we 
have that |5(i?BEc)| — vV^- Then we can evaluate with 
the help of the formula in Eq. H135|l the binding energy 
for this value of the magnetic field. This gives 



em(i3BEc) — -vV^ 




(136) 



Substituting this expression in the definition of Z in 
Eq. (|52() and expanding the result in the small param- 
eter rj^ /tp, we obtain 



^(-Sbec) — 1 



1/4 



1. 



(137) 



Therefore, for a narrow resonance the bare component 
of the dressed molecule cannot be neglected. It turns 
out that also in the so-called universal region, where 
l/kF\a{B)\ ^ 1, we have {) < Z < 1 and the thermo- 
dynamics is actually not universal. 

For a very broad resonance, the bare component Z 
in the universal region l/kF\a{B)\ ^ 1, is not strictly 
zero but it is rather small. When the chemical potential 
approaches half the binding energy, the bare component 
weight is still Z ~ 0. Infact, expanding Eq. H136|l in 
tpjr^ we find 



ZiBBEc) ^ 24 < 1. 



(138) 



This implies that at magnetic fields where the BEC-BCS 
crossover takes place, the dressed molecule is almost com- 
pletely in the open channel component. Furthermore, the 
binding energy in Eq. H135|) can be approximated like in 
Eq. ^ by 



ma'^{B) 



cx 6" 



(139) 



which is equivalent to neglecting the term 2/i in Eq. (|134|l . 
Note that this is the same solution determined by the 
pole of the the many-body T-matrix in the single- channel 
approximation of Eq. (|77|l . It is also clear now why, in the 
case of a broad resonance, it is essential to replace in the 
equation of state the detuning with the curve of Eq. H127() 
in order to avoid inconsistency in the BEC limit between 
the two crossover equations. In conclusion, in the case of 
a broad resonance, we can roughly speaking distinguish 
between two different crossovers. At detunings where 
the many-body crossover to a Bose-Einstein condensate 
of dressed molecules takes place, the dressed molecules 
live essentially in the states of the atomic continuum. At 



larger negative detuning, the wavefunction of the dressed 
molecules undergoes a second two-body crossover to the 
bare closed channel eigenstate. The energy scale of this 
second crossover is 77^ as mentioned above. 

In general, as far as the background interaction satis- 
fies the weak-coupling condition feFl^bgl < 1, the above 
discussion does not change qualitatively, if we extend the 
analysis of the many-body T-matrix in Eq. (|123|l without 
taking the limit Obg = 0. The gap equation in Eq. H130() 
in the BCS weak coupling regime modifies to 
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Tg^(O) 



1 



2^ - SiB) 



S 0,0^ 



0. (140) 



In the case of a broad resonance the weak-coupling solu- 
tion is again given by Eq. H132f) and can now be written 
as 
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Tc ~ -r-^e'^ ^ exp 



2kF (kr- 
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(141) 



At large negative detuning the pole of the many-body 
T-matrix (|123|l reproduces the correct two-body physics 
of the dressed molecules. The equation for the chemical 
potential in Eq. H134() assumes the general form as in 
Eq. gni) 



5{B) - 2/i 
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-2fj,7n 
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Therefore the chemical potential goes to half the bind- 
ing energy of the dressed molecules em given by the more 
general solution (|47|l . The analysis concerning the dif- 
ferencies between the broad and narrow resonances re- 
mains unchanged, except that the energy of the second 
two-body crossover defined above is now near /ma^^. 
This analysis holds in general for the resonances currently 
used in experiments. One exception, however, concerns 
the region at large positive detuning of the the broad res- 
onance at 834 G in a gas of of ^Li atoms. In that case 
the background interactions turn out to be so large that 
h^/ma^g ~ ep- In other words, the gas never enters a 
weak-coupling regime at positive detuning and the gap 
equation in Eq. H140() is not expected to be valid. 

Figure [201 shows the mean-field curve for the critical 
temperature for a gas of ""^K atoms near the broad res- 
onance near 201.2 G. The mean-field approach based on 
Eqs. (|125|l and (|128|l does not describe a smooth BEC- 
BCS crossover because it neglects finite lifetime effects 
of the molecules. The curve for Tc exibits a sharp phase 
transition between a Bose-Einstein condensation phase of 
dressed molecules with binding energy Cm to a Bardeen- 
Cooper-Schriffer superfluid. This phase transition is only 
an artifact of the mean-field aproximation because the 
two phases break the same U{\)/Z2 symmetry |ll] . How- 
ever, the dressed-molecules are the real physical compo- 
nents of the system. Therefore, the location of the kink 
in curve for Tc of Fig. 1201 assumes an interesting inter- 
pretation. According to the poor man's approach of Ref. 
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Figure 20: Phase diagram of atomic ''"K near the broad res- 
onance near 201.2 G, as a function of magnetic field and tem- 
perature for a Fermi temperature of the gas of Tp = 0.35 jjK. 
The solid line exhibiting a kink at about 0.5 G gives the crit- 
ical temperature for either a Bose-Einstein condensation of 
dressed molecules or a Bose-Einstein condensation of atomic 
Cooper pairs. The critical temperature for the latter is cal- 
culated by simultaneously solving the BCS gap equation of 
Eq. 112511 and the equation of state of an ideal mixture of 
atoms and dressed molecules as in Eq. lEHl- For compari- 
son, the upper dashed curve is the analytical BCS result of 
Eq. I132II . The lower dashed line is the crossover between the 
two Bose-Einstein condensed phases. The other solid smooth 
line describes the result of the Nozieres-Schmitt Rink approx- 
imation. 



[43, it fixes a scale in the magnetic field below which, 
the superfluid transition temperature is determined by 
the condensation temperature of the dressed molecules. 
Note that this value of the magnetic field coincides also 
with the magnetic field value at which the size of the 
dressed-molecule, that is equal to a{B), approaches the 
interparticle separation kp^ 37]. It is important to stress 
that, without the replacement in the equation of state of 
the detuning with the true location of the molecular level 
in presence of coupling described by Eq. (|127|l , this scale 
would have been orders of magnitudes smaller. Interest- 
ingly, the location of the crossover line in our mean-field 
picture that is shown in the dashed curve of Fig. |2D1 co- 
incides with the threshold for the non zero condensate 
fraction in the (T/Tc, AB) plane found by Regal et al. 
in the experiment of Ref.|l5j| that we have discussed in 
Sec.lnl 



C. Theory of Gaussian fluctuations. 

A smooth curve for the critical temperature Tc can 
be obtained by including finite lifetime effects of the 
molecular field in the theory. This can be achieved by 
considering the quadratic fluctuations on top of mean- 
field, according to the Nozieres-Schmitt Rink approach 
[l2,llal45ll47ll50| . In the Nozieres-Schmitt Rink approx- 



imation, one equation is still determined by the Thouless 
condition expressed by Eq. (|140|l . The equation of state, 
however, is modified because it has to incorporate the 
effects of Gaussian fluctuations around the saddle point. 
The new equation of state reads ,47j| 

+ ^Tr[l„(;3G-')]), 

where the full molecular propagator G is given by 
Eqs. ijHll) and ijHSI) and the symbol Tr[ ] denotes the 
sum over the Matsubara frequencies and momenta. This 
equation describes a gas of atoms and dressed molecules 
[ll8j . Using the definition of the full molecular propaga- 
tor Eq. H143() can be rewritten as 

n = (144) 
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The partition of the different contributions in this equa- 
tion of state suggests a possible interpretation of each 
term on the right-hand-side of Eq. 11441) . The first term 
denotes the number of unbound atoms. The second term 
describes the contribution of the bare molecular compo- 
nent while the third one defines the open-channel com- 
ponent of the dressed molecules. The analysis of the 
BEC extreme of the crossover reveals that this interpre- 
tation is correct in that regime. In the BEC limit the 
chemical potential approaches half the binding energy of 
the dressed molecules, i.e., /i ~ em/2. Therefore, the 
chemical potential becomes large and negative and the 
contribution of the unbound atoms can be neglected. In 
first approximation the gas constituted mainly of dressed 
molecules 



71 ~ 2 

We want to show that in that limit 
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Taking the contribution of the pole of the dressed molec- 
ular propagator we have 
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where we have used the definition of the residue Z 
given in Eq. (|52|) . and that in the BEC hmit we have 
i9/9/i ~ 2d/dE because the many-body molecular self- 
energy approaches its two-body limit given by Eq. H43|) . 

In the remaining part of this section we want to show 
that in the limit of a very broad resonance our two- 
channel model reproduces the results for Tc found by Sa 
De Melo et al. in Ref. based on a single- channel 

model nT 



As we have shown in Sec. lIII El in the case of 

a broad resonance, the general expression for the pairing 
field propagator /iGa given in Eq. H71(l reduces to that 
given in Eq. H78|) . The latter corresponds exactly to that 
defined by Sa De Melo et al. in Ref. 0| . This argument 
is sufficient to proof that our two-channel model is able 
to reproduce the results of the single-channel model in 
the limit of a very broad resonance when if' 3> cf ■ How- 
ever, a detailed analysis of our crossover equations for 
a broad resonance, gives us more insight about the two 
different descriptions. When if ^ ep, the gap equation 
in Eq. (|140|l can be approximated with 



1+ (T2B(o)+5 
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S(O,0+)~0, (148) 



along the all range of magnetic fields where the BCS-BEC 
crossover takes place. This is possible because, as we have 
explain in Sec. IIII El in the case of a broad resonance 
we have \S{B)\ \2fi\,em except of course very close 
to resonance. However, in that region the gas becomes 
essentially universal and it can be shown that the non 
universal correction term 2/i can also be neglected 38] . 
By using Eq. iB^ Eq. 1] 148(1 can be rewritten as 



i + (^^-^js(o,o-).o, 



(149) 



which corresponds to the single-channel gap equation of 
Sa De Melo et al. in Ref. 13]. A similar analysis con- 
cerns the equation of state in Eq. H144() . In the case 
of a broad resonance the weight of the bare molecular 
component remains quite small such that Z ~ even 
when the chemical potential approaches half the binding 
energy ~ Em — f? jma^ . Therefore we can neglect com- 
pletely the bare molecular contribution to the number of 
particles and the equation of state reduces to 
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(150) 



Furthermore, in the crossover near a broad resonance the 
bare dynamic of the molecular field can be neglected and 
the molecular propagator in Eq. (|150l) can be approxi- 
mated as in Eq. H75|l . Therefore Eq. 1)150(1 can be rewitten 
as 



1 



1 



'h(3V 



ETr[G.„]+^Tr 



d 



-5{B) - ^^l 



(151) 



d_ 



(2^ETr[ln(G-;,)]- 



(152) 



(3V 



Tr 



In ( 



-5{B) - SE„ 



At this point, to make this part of the derivation more 
transparent we consider the special case when the Obg — > 
0. In that case the molecular selfenergy of Eq. I(83() is 



fiEm (q,w„) =5^-^E 
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Using the expression of the total scattering length in 
Eq. 144() this can be rewritten as in Ref. ^3 

" = |(^ETr[ln(0]+^Tr[ln(G.-]), 

(154) 

where Ga is the pairing field propagator in the Gaus- 
sian approximation defined in Eq. ((78(1 . In the general 
case when Obg is finite, the derivation is analogous but 
only more involved because the energy dependence of the 
atom-molecule coupling 119]. 

By comparing Eq. 1(150(1 and Eq. ((154(1 we argue that 
the propagator Ga of the composite boson in the single- 
channel picture contains the open-channel contribution 
of the dressed-molecules. The precise mathematical rela- 
tion between the two entities will be discussed further in 
one of the following sections where we a detailed analysis 
of the BEC limit in the superfluid state. 

The curve of the critical temperature Tc for the broad 
resonance in atoms based on the crossover Eqs. ((149(1 
and l(154() is shown also in Fig. [201 This result shows 
that our modified mean-field approach, based on the two- 
body physics of the dressed molecules, can be considered 
a rather good approximation in view of its simplicity. 



D. Self-consistent equations. 

The curve of Tc in Fig. [201 calculated according the 
Nozieres-Sclimitt Rink approximation exhibits a maxi- 
mum. This maximum disappears in a self consistent 
approach The latter includes some fluctuations 

that lead to an increase of the effective mass of the 
dressed molecules and therefore to a descrease of the 
transition temperature. A self-consistent approach for 
the calculation of Tc in a dilute gas of fermions inter- 
acting with strong attractive interactions has been de- 
veloped by Haussmann in Ref. ^14] in relation with a 
single-channel model and has been extended by Xiang 
and Ji in Ref. to the two-channel Feshbach model. 
Their numerical results shows clearly that the curve for 
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Tc in the BEC-BCS crossover evolves monotonically from 
one limit to the other when selfenergy effects between 
atoms and molecules are treated self-consistently. In a 
self-consistent approach, the conservation of the number 
of total particle in the BEC-BCS crossover is given by 
the equation 



Therefore the equation of state in Eq. (|155|l can be rewrit- 
ten as 
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where the full atomic Green's function is calculated from 
the many-body T— matrix as in Eq. H109|) but including 
to all order the selfenergy effects in the fcrmionic and 
molecular propagators. In this section we want to discuss 
the two coupled Eqs. (|155|l and H140|) in connection with 
our dressed-molecules picture. The analysis of the BEC 
limit is particularly illuminating to this purpose. We 
consider here the case when the background interaction 
can be neglected. However, as we have seen, it is always 
very easy to generalize our result including the effect of 
the non resonant interactions. In the BEC limit we have 
2/1 ~ Em and ^ ksT^ep and the fermionic selfenergy 

can be approximated to the lowest order in the fermionic 
propagator as Q 

ni]^(k,c.„) ^-G_,,o(k,c.„) 52_i_y'rf3k'_L^G(k',c.;), 

(157) 

Z n^{T), (158) 
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Far enough in the BEC limit, it is always possible to 
expand the first term on the right hand side in powers of 
^Pg/I^l - Apg/em and we get 
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where Z in the BEC limit is g iven by taking the limit 
Obg ^ in Eq. ifS^ . that is 



Z = 



2^/k„ 



(159) 



This means that the fermionic propagator assumes a dou- 
ble fraction structure as 



that can be rewritten as 

n = 2 (1 - Z) n^{T) + 2 Z n^{T) = 2 n^{T). (164) 

Therefore, at sufficiently negative detuning we find a gas 
of thermal dressed molecules. However, we have not spec- 
ified the approximation of the molecular selfenergy in 
the molecular propagator G. To the lowest order, the 
molecular selfenergy HT,,^ is given by the expression in 
Eq. H83|) calculated taking two free fermionic propaga- 
tors Gcr,oG-a;o- In this approximation, Eq. (|164|l repro- 
duces the gas of non interacting dressed molecules as in 
the Nozieres-Schmitt Rink approach found in Eq. (|146|l . 
The interaction between the dressed molecules at T > Tc 
is introduced calculating the molecular selfenergy to the 
lowest order in the fermionic selfenergy. This requires 
to take one free and one dressed fermionic propaga- 
tors (Gcr,oG_cr) in the diagram of the molecular selfen- 
ergy. Taking, in the lowest order, the selfenergy in the 
fermionic propagator given by Eq. (|156|) 

/^I^^(k,^„)c.-(^/d3k'^E (165) 

G(k + k',c^„+c^:j G_„,o(k',c^;), 
the molecular propagator assumes the form 
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that is equivalent to a fermionic BCS propagator charac- 
terized by a "gap" in the spectrum defined as 
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The first-order correction hAT,^ is given by 
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with the many-body molecule-molecule vertex given by 
r^°™ (K + K', K + K', (K + K')/2, {n„ + i^n)/2) = 



- w„). 
(168) 



In the BEC limit at temperature near Tc, we have 
hn{B)\ > Ml.h^K^lAm and lAE^I < |Sm| and the 
molecular propagator of Eq H166() can be expanded 
around the pole of the molecular propagator without the 
/lASiji correction. This gives 
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where we have defined the dressed-molecule chemical po- 
tential as /ini = 2/i — e-ax{B). Furthermore, the zero- 
momentum and energy correction to the molecular self- 
energy in Eq. p69|) is given by 
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Evaluating the Matsubara summation and the momen- 
tum integral of Eq. (|168|1 in the BEC limit where 2fj, ~ em 
we find 
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where we have used the definition of Z in Eq. (|159|l . 
In the case of a very broad resonance at large nega- 
tive detuning, when the gas of dressed molecules enters 
the weak-coupling regime kFa{B) < 1 and the chemi- 
cal potential is essentially equal to half the binding en- 
ergy of the two-body dressed molecules, the wavefunc- 
tion in Eq. (|53|l of the dressed molecule can still contain 
only a very small amplitude in the closed bare molecular 
channel. This means that the binding energy is in the 
quadratic regime in the magnetic field given by Wigner's 
formula Cm — —h^/ma^{B) and that Z ~ 0. In this case 
Eq. I|172l) can be written as 

Z /lASm (0, 0) ~ 2 ^''^'^"^'^^ n^T) (173) 
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This term in the dressed molecular propagator of 
Eq. I|169|l represents a many-body mean-field shift of a di- 
lute gas of molecules interacting with positive scattering 
length 2a{B). 

The interaction energy of the dressed molecular gas 
has been measured 18j by fitting a Thomas-Fermi pro- 
files to the spatial distributions of the trapped ultra- 
cold Fermi gas at densities of about 10^^ cm"^ in the 
BEC limit. These experiments consider the crossover 
across the broad resonance at about 834 G in atomic 
^Li atoms. They assume, for magnetic fields between 
600 G and 780 G, a dilute gas of tightly-bound inter- 
acting molecules. They find that the molecule-molecule 
scattering length is in agreement with am{B) ~ 0.6 a{B). 
The origin of the value of 0.6 a{B) has been explained 
by solving directly the Schroedinger equation of the four- 
body problem 1761 and also in terms of Feynmann dia- 
grams in Refs. |39l l40l | by using a single-channel model. 
In the latter approach, the authors of Refs. J9, 40] have 
identified the exact equation for the generalized molecule- 
molecule vertex that contains the dominant fluctu- 
ations in the BEC limit leading to the value found in 
the experiment. Considering the same class of diagrams 
in our two-channel formulation for the vertex Fm intro- 
duced in Eq. H168|l would also lead to a dressed molecule- 
molecule T-matrix F^ (0,0, 0,0) ~ iTTfi^ 0.6 a{B)/2m. 
For example, summing over all the ladder diagrams in 
the particle-particle channel beyond the Born approxi- 
mation described by Fj^"'''^ one would recover the result 
F,„ (0,0,0,0) ~ Airh^ 0.75 a{B)/2m found by Fieri and 
Strinati in ;33|- 

The interaction between dressed molecules introduced 
in Eq. Ij 173(1 does not cause any shift on the critical tem- 
perature of the Bose-Einstein condensation of the dressed 
molecules because it does not shift the efi'ective mass of 
the dressed molecules. Therefore, we have to consider 
the momentum and frequency dependence in the inter- 
action vertex F^"'" in Eq. ((168(1 . Follo wing the analysis 
of the BEC limit given by Haussmann [lj]j we consider 
the ansatz 
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(174) 



for the drcssed-molecular propagator. The Nozieres- 
Schmitt Rink approximation corresponds clearly to the 
case Zi = 1 and m* = m. To find the first order correc- 
tions to Zi and m* , we substitute the ansatz of Eq. 1(174(1 
in Eqs. ((1671169(1 . we perform the frequency sum and the 
momentum integral in Eq. ((168(1 and finally we expand 
the resulting expression up to the first order in and 
f2„. Ultimately, in the case of a broad resonance, we find 



m 
Zi 

Mm 



m [l + 2n {n,na{B)^)] (175) 
1 + 6tt n^a{B)^ (176) 
2FB-"(0,0,0,0)n,n(r). (177) 



The critical temperature in the BEC limit is determined 
by the equation of state. In a broad resonance, one gen- 
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erally has still Z ~ when the gas enters the BEC limit 
of the crossover. In that case, the number of particles in 
the equation of state of Eq. (| 15 5(1 is determined mainly 
by the fermionic Green's function 
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(178) 



The evaluation of this term given in Eq. ((164|l is calcu- 
lated taking the fermionic propagator only the first- 
order in the selfenergy. Expanding the fermionic Green's 
function at the second order in its selfenergy, one can 
show 14] that the result of Eq. (|164|l is modified as 
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Substituting in this equation the propagator of Eq. 11741) 
the Z and Zi factor drop out. Therefore the equa- 
tion of state describes effectively an ideal gas of dressed 
molecules with renormalized mass m* . The critical tem- 
perature is calculated from the relation 
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This leads to the negative shift of the critical temperature 
Tc respect to the ideal Bose-Einstein critical temperature 
Tbec 

ATc/TBEC^m/m*-l^- ^ (181) 

(3tt) (kpa(lJ)) 

found by Haussmann. 

Although the self consistent calculation is more gen- 
eral than the Nozieres-Schmitt Rink approximation, we 
think that the correct theory should lead to a positive 
shift and then to a maximum in the curve of Tc- Infact, 
accor ding t o the theo ry of dilute weakly-interacting Bose 
gases |l20l Il2ll Il22 ^ , the repulsive interactions between 
the molecules in asymptotic BEC limit are expected to 
enhance the value of the critical temperature. The phys- 
ical origin of this maximum, therefore, is different than 
that described by the Nozieres-Schmitt Rink approach, 
as the latter neglects the interactions between the non- 
condensed molecules. 



V. BEC-BCS CROSSOVER BELOW Tc. 

In this section, we consider the BEC-BCS crossover in 
the broken symmetry state. The Gross-Pitaevskii theory 
at T = of the dressed-molecule condensate and the Bo- 
goliubov theory at finite temperature are derived in the 
BEC limit by means of analytical methods. As before, 
the single channel results are derived as a special case in 
the limit of very broad resonances. Moreover, the relation 
between the composite boson of the single-channel and 
the dressed-molecule of the two-channel model is clearly 
formulated. 



A. Mean-field theory at T = 0. 

At T = a set of two self-consistent equations for the 
BCS energy gap and the chemical potential can be de- 
rived generalizing the variational approach introduced by 
Leggett ^3 to our two-channel atom-molecule Hamilto- 
nian given by Eq. H2UI) . The zero-temperature ground 
state can be written ^ ,55,] as a product state of both a 
bare molecular and an atomic contribution, i.e.. 
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where the normalized fermionic wavefunction is the stan- 
dard crossover ground state 0, ^ 



K> = nK + ^;kcIcl,)|0), 
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and the normalized bare molecular contribution 1^"^) is 
given by 



(184) 



The real variational parameters are Uk, I'k and y/riQ. 
Minimizing the grand-canonical energy with respect to 
the parameter i/uq, we find the constraint 
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where Abg is the order parameter related to the attractive 
background interaction 

Abg = XI VbgUkWk = X! Vbg(a-kiak^i) ■ 

k k 

The variation with respect to v^, under the normalization 
constraint f k + "k = 1 > yields 
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'(ek-M)^ + |Ap 
where the total gap of the theory A is defined as 

A = Abg - .gbarcV^- (188) 

Substituting these results in the definition in Eq. H186|l 
and using the constraint in Eq. 1(185(1 . we arrive at the 
gap equation in the compact form 
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Figure 21: BCS energy gap and chemical potential as a function of the magnetic field for the very broad resonance at 834 G 
in a ®Li mixture at low density. 




Figure 22: BCS energy gap and chemical potential as a function of the magnetic field for the very broad resonance at 834 G 
in a ®Li mixture at high density. 



This equation can be rewritten in terms of the renornial- 
ized quantities as 
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as a consequence of the general result for the many-body 
T— matrix in Eq. H33|l . combined with the fact that the 
gap equation in Eq. p9()|l must be equivalent to a sum 
over all ladder diagrams. The description of the BEC- 
BCS crossover at the mean-field level of approximation 
is completed by the equation of state 
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The Eqs. ((M]|l - lfTI?^ and ifTn^)) represent a set of cou- 
pled equations in A and /i at fixed density. Note that in 
contrast to the normal state result, the analysis at zero 
temperature in the broken symmetry state does not re- 
quire any modification of the equation of state in order to 
get the proper two-body physics of the dressed molecules 
in the BEC limit. This is because the equation of state 
of BCS mean-field theory at T = sums automatically 
over the ladder diagrams. As a result, we expect also to 
have smooth crossover curves for A and fi already at the 
mean-field level. 

At positive detuning, in the weak-coupling regime 
kF\a{B)\ < 1, the binding is a cooperative effect in the 
vicinity of the Fermi surface. The gas consists of largely 
overlapping weakly-bound Cooper pairs. The chemical 
potential is fixed at /i = e^? by the equation of state in 
Eq. H193() . The solution of Eq. H190() is then easy calcu- 
lated analytically with the result 



Sep _2 2 
A ~ - — e 6" 1 exp 



2kF\a{B) 



(194) 



The energy gap is related to the BCS transition critical 
temperature H131|) by the relation Tc — [e? I'k) A, as is 
expected for weak-coupling superconductivity. Note that 
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the energy gap of Eq. (|194|l coincides formally with the 
shift in the ground state energy of Eq. H97|) . This is not 
surprising because the singularity in molecular selfenergy 
discussed in section IIII Fl and the Cooper pairing occur 
in the same diagram. 

At negative detuning, the system evolves to a Bose- 
Einstein condensate of tightly bound dressed molecules. 
The roles of the gap equation and the number equation 
are exchanged. In the BEC limit, the chemical potential 
becomes large and negative. The crossover equations can 
be expanded in powers of |A|/|/i| as 



T 



2B 



5-2/1 



2m VIMI 



16 /i2 



n =2no + |A| 



and 
(195) 
(196) 



Neglecting the term quadratic in |A|/|/i|, Eq. p95|l be- 
comes equivalent to Eq. H142I) . Thus the chemical poten- 
tial approaches half the binding energy of Eq. H47|l of a 
Bose-Einstein condensed dressed molecule, as expected. 

The curves for fi and A obtained by solving numeri- 
cally the mean-field crossover Eqs. (|19U|l - (|192f) and H193|) 
are shown in Figs. [2] and [53 for a gas of ^Li atoms at 
two different densities and in Fig. [221 for a gas of 
atoms. In both cases we have ry^ S> ef- In each figure 
the vertical lines represent the boundary of the strong 
couphng regime kp\a{B)\ > 1. In the ^Li gas with high 
density of Fig. [22 the gas never enters the weak-coupling 
regime kp\a{B)\ < 1 in the range of the magnetic field 
considered in the figure. 

The mean-field curves for A and /i do not differ quanti- 
tatively from the solutions obtained by the single-channel 



approximation pTl [T^ based on the equations 
1 1 



rp2B _ ii iirh^ajB) 
bg 5 m 



(197) 



-y 

V ^ 



^ k V2V(ek-A^)' + |A3cP) 



2ek 



and 



(27r)3 



fk - Ai 



(ek - 



(198) 



where the single-channel gap is defined as Asc = (Vbg + 
Vrcs)(V-'T(^)''/'i(^))- These equations are obtained from 
Eqs. H190() and (|193|l neglecting explicitly the contribu- 
tion of the bare molecular boson. They can be derived 
also by integrating out the molecular field in the action 
of Eq. H60|l from the very beginning. 

At negative detuning the chemical potential goes to 
half the binding energy ~ —h^/raa?{B) described in 
Figs[2IIb),[22Ib) and[^b) by dotted lines. In Fig.[2Si;b) 
the curve of the binding energy of the dressed molecules 
is continued at positive detuning by the location of the 
maximum in the two-body density of states given in 
Eq. I|58|) . This illustrates the main idea of our modi- 
fied mean-field picture of Sec. IIVBI The magnetic field 
at which this curve crosses 2^ is indicated by a vertical 
line. At this magnetic field the energy of the two-body 
dressed molecules approaches the chemical potential of 
the gas [ill. 

In Figs. [^ [221 and [23 the dotted lines represent the 
two extreme limits of the crossover. The lines denoted 
as BCS correspond to the BCS solution of Eq. (|190|l . 
T he lines denot ed as BE C correspond to the curve A = 
■\/ (16/37r)e p / ykpaiB), which is the asymptotic solution 
for the energy gap according to the single-channel gap 
equation in Eq. H197|) . 
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B. Asymptotic limit of Z in the deep BEC limit. 

The analysis of the equation of state in the BEC Umit 
requires some special consideration. The correct picture 
of the crossover across the Feshbach resonance is based on 
the dressed molecule, which is the true energy eigenstate 
of the diatomic molecule in the presence of the atom- 
molecule coupling. Unfortunately, however, the dressed 
molecule density does not appear explicitly in the cou- 
pled equations of our mean-field treatment and mean- 
field theory is therefore not able to extract this quantity. 

Nevertheless, a gas of weakly-interacting closed- 
channel bare molecules is realized at sufiiciently large 
negative detunings. More precisely, the total density 
amounts to 
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n ~ 2no, (199) 

when kpa{B) < 1 and > h'^/ma'^^. For less nega- 

tive detunings we must use ;54i] 

no = Zn^c, (200) 

where Time is defined as the condensate density of the 
dressed molecules. A careful analysis of the equation of 
state in Eq. 119t)|) shows that this definition is indeed 
consistent. The proof proceeds as follows. First we as- 
sume that in the deep BEC limit the total density can 
be approximated by the density of the dressed molecules 

n ~ 2n„,c, (201) 

with binding energy Cm — 2/z given by Eq. (|T7|) . Then, 
we eliminate the explicit dependence on S and Abg by 
using Eqs. H192fl and H142|l . to rewrite the equation of 
state as 

2n^, ~ 2Zn,„c + ^^I^ _ (202) 

Solving for Z we recover the correct two-body limit 
expression of Eq. (|52|l . i.e., 

r n -1 




Putting back this Z in the equation of state we find again 
Eq. P7T|) . 

It is important to stress the fact that the approxima- 
tion introduced in Eq. (|201|l is valid only deep in the 
BEC limit of strongly-bound dressed molecules. This 
can be understood as follows. According to the BCS 
theory, at zero temperature all the atoms are paired 
far below resonance. In our case we have a gas of 
dressed molecules which share an amplitude on the open- 
channel Cooper pairs wave-function and on the bare 
closed-channel molecular bound state. However, even at 



Figure 24: The value of Z measured in j20( | is plotted as func- 
tion of the magnetic field (dots). The continous line repre- 
sents the two-body Z given by Eq. H52^ The vertical dashed 
lines represent the boundary of the strong coupling regime 
kF\a{B)\ > 1. Below 800 Gauss, we have 2/i ~ em and the 
two-body Z matches perfectly with the experimental value. 
The Fermi energy of the gas is at Tf — SSOnK. 

zero temperature not all the dressed molecules are con- 
densed because interactions deplete the condensate den- 
sity as the interactions between the Cooper pairs do in 
the zero-temperature BCS theory. Therefore, Eq. H201|) 
neglects the depletion of the dressed-molecules conden- 
sate and is correct only asymptotically in the deep BEC 
limit when the depletion descreases progressively because 
the gas becomes more and more dilute. 

This means that Eq. (|201|l is valid only at magnetic 
fields such that kpa(B) < 1, when the many-body 
part of the BEC-BCS crossover has already taken place, 
and where the chemical potential approaches half the 
molecular two-body binding energy in Eq. H47|l . When 
kF\a{B)\ > 1, the fluctuations in the atomic fields be- 
comes important and also at T = the gas cannot be 
approximated only by a condensate of dressed molecules 
|54l |. As a result, the bare molecular probability Z is 
not given by 2no/n [5^ Il23| | as in the asymptotic BEC 
limit. Moreover, Z deviates from its two-body expression 
of Eq. (|52() and the two mean-field equations are not suffi- 
cient to calculate its value, because they can not calculate 
the depletion of condensate of the dressed-molecules. A 
more refine approach is needed beyond the saddle point 
solution [s^ . The quantity Z has been recently mea- 
sured by Partridge et al. Their data are illustrated 
in Fig. 1241 The experimental data are in very good agree- 
ment with our two-body Z of Eq. H52f) in the BEC limit. 
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C. Gross-Pitaevskii theory of the dressed-molecule 
condensate at T = 0. 

We show in this section that, in contrast to the nor- 
mal state analysis, the Born interactions between the 
dressed molecules in the BEC limit of a broad resonance 
described in the previous section, arises already at the 
mean-field level in the crossover theory below Tc based 
on Eqs. H19U|) and H193|) . The effective chemical potential 
for a dilute weakly-interacting gas of dressed molecules, 
is defined, at the leading order, as 



(203) 



where am is scattering length of two dressed molecules 
with mass mm — 2m. The quantity /im is given by 



2^ 



(204) 



where 2/i is calculated retaining the quadratic term 
|Ap/|/ip in the gap equation in Eq. H195|) . Using 
Eqs. H200() . (|201|l and that 2/i ~ em the gap equation 
in Eq. p95|l can be rewritten as 



47r?iVff(B) 
2/i -I- rime ^- 



or 



fJ>Yn — ^mc ^ 



m 



(205) 



r, (206) 



where the binding energy em is given by the general solu- 
tion of Eq. (|47|l , and the quantity Qcs has the dimensions 
of a length and is given by 



aeff(-B) 



(207) 



Using in Eq. (|2U5|I the definition of Z in Eq. we 
find for the effective chemical potential of the dressed 
molecules at the mean-field level 



Mm ^ "-mc 2acff(B)(l - Z). 

mm 



that implies 



2a,s{B)il-Z). 



(208) 



(209) 



However, for very broad resonances such as that at 834G 
in the ^Li gas, when the gas enters the weak-coupling 
regime kpa{B) < 1, the wavefunction in Eq. (|53|l of the 
dressed molecule still contains only a very small ampli- 
tude in the closed bare molecular channel. This implies 
that the dependence of the binding energy on the mag- 
netic field is still in the quadratic regime by Wigner's 



formula em — —h^/ma'^{B). Therefore, in this regime 
Eq. H208I) can be approximated as 



nmc^ 2a{B){l-Z). 



(210) 



This re sult is eq uivalent to that of a Gross-Pitaevskii 
theory |l24l Il25l | for a gas of Bose-Einstein condensed 
dressed molecules interacting with an effective scattering 
length 



2a{B)(l - Z). 



(211) 



However, because in this regime we have Z ~ 0, the 
scattering length between the dressed molecules in the 
Born approximation in the mean-field approximation is 
given essentially by 

a^°™ ~ 2a{B). (212) 
Such a situation is typical of a very broad resonance. 

D. Dressed molecules versus composite bosons. 

Our derivation in the previous section, has to be con- 
sidered as an extension of the mean-field analysis of the 
BEC limit for a single-channel model crossover discussed 
by Sa de Melo et al. in Ref. [3 (see also Ref. jH] ) based 
on Eqs. H197|l and H198|l . In the mean-field single-channel 
model the Cooper pairs evolve in the BEC limit towards 
a dilute gas of weakly-interacting composite bosons with 
scattering length 2a{B). In the BEC limit, the equation 
of state of the mean-field single-channel model at zero 
temperature can be approximated by |l4l IsM l35l | 



m^a{B) 

471^4 I 



(213) 



The composite boson condensate of the single- 
channel model is defined from [l4L l33L ISal 



/ 87rft4 
'm?a{B) 



(214) 



such that the BEC limit of Eq. I|213(l is characterized by 
a gas of composite bosons 



21* 



B\2 



(215) 



The equation of state in Eq. (|213|l must be compared 
with its two-channel analogue in the BEC regime where 
2/i ~ Cm- We have 



n ~ 2no -I- |A 

= 2Z7imc 



ATTy/2^^\h? 

2Zn,^c'r] 



(216) 




= 2n„ 
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Moreover, for a very broad resonance, we have 2/x ~ em — 
— h^/ma^{B) and Z ~ 0. Therefore it is ^ Cbg, and 
we find 



3 



(217) 



This relation between the residues of the poles asso- 
ciated with the composite boson and with the dressed 
molecule of the two-channel model in the case of a 
very broad resonance can now be better understood if 
we remember the phy sical meaning of Z |46l | and of 
Snfi^/m^aiB) [ll (ll |3|. The factor 8nh'^/m^a{B) 
rcficcts the difference between the bosonic propagator 
of the single-channel composite boson and the particle- 
particle ladder propagator in Eq. H78|) in the BEC limit. 
In that limit we have 2\^i\ ~ \h^/ma{B)^\ > feT and 
the latter can be approximated as 



(218) 



a{B) 




4m 



2^ 



Expanding for small a{B) this can be rewritten in the 
polar form 35] 

Ga (K, f]„) /;i=T^B (j^^ (219) 

87rfi4 1 



T?a{B) ihnr. 



4m 



MB 



where the chemical potential of the composite boson is 
defined as fiB = 2^ — H^/ma{B)'^ and can be calculated 
from the single-channel gap equation H197|l in the BEC 
limit. We find [11 1^3 



Ms 



ATTh^2a{B) |As 



2to 



4|/i| 



(220) 



To be rigorous, we should have considered the particle- 
particle ladder propagator in the superfluid state. How- 
ever, in the BEC limit, hftn, ?i^K/4m and mb are of the 
same order or smaller than h^/ma{B)^. Therefore the 
mean- field inverse particle-particle propagator hG~^ , up 
to the first order in hfln and fi^H/Am, is the same as in 
the normal state. 

The factor Z refiects the difference between the bare 
and the dressed molecule |4^J . The full molecular propa- 
gator in Eq. 

hG^^ (K, Vtn) = ihQn - -. — + 2^ - 6iB)+ 



1 



-> |g-(0,2ek) 



V ^ 

k 



4rn 
1 - 



K./2+'k 



K/2-k 



ihflr. 



2n - 2ek - eK/2 



1 
2^ 
(221 

can be approximated in the BEC regime, where 2/i 2 
£„! — —h^/raa?{B) and Z ~ 0, as 



hG 
1 



ihQr 



4m 



-> |g-(0,2ek) 



V ^ 

k 



2^1 ~ 6(B)- 
1 



2/i - 2ek - eK/2 



1 
2^ 

(222) 



However, for broad resonances, in this limit we have 
also that \S{B)\ » 2|/i|, fi^K/Am and hn^, fi^K/Am < 
hi] (K, il„). Therefore, we can approximate the dressed- 
molecular propagator further as 



hG-^{K,nn) ~ -6iB) 
-l5:i5^^0,2eOP 



1 



(223) 
_^ 1 

'2m — 2ek — eK/2 2ek 
(224) 



Calculating the integral this can be rewritten as 



(225) 



a{B) 



^^j-ih^n 



h^K 

4m 



- 2/i 



from which, expanding as in Eq. 



we obtain 



ft-iG(K,ao- ^ 



&Trh'^ 



g2 m^a{B) ihn^-!^ 



(226) 



where Mm = Mb- However, expanding Eq. (|223|l around 
the real pole at Cm — h^ /ma^{B) the propagator can also 
be written as 



h-'G{K,n^)~z{B) - 



(227) 



Comparing Eqs. (|226|) and H227|l we find again the re- 
lation in Eq. H217|) . Our above analysis is based on a 
mean-field approach and due to this reason the propaga- 
tors of Eq. (|226|) and H219|) describe free bosons. This is 
because we are considering only the bosonic propagators 
only up to first order in hfln and h^Ji/Am. 



E. Gaussian fluctuations around the saddle-point 
solution. 



In order to get in the BEC limit a superfluid of in- 
teracting bosons, we have to extend our analysis beyond 
the saddle point approximation. The nonzero tempera- 
ture molecular propagator in the broken-symmetry state 
at the level of Gaussian fluctuations can be approximated 
in the BEC limit by the matrix 
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_ 1 ( -ifiujn-2^i + e^/2 + 5{B) \ / ^En (K, cj„) ^Eia (K, cj„) \ 

' zfii^„-2/i + eK/2 + ^(B) i + UE2i(K,c^„) nE22(K,c^„) ^ ^^^^^ 



where the molecular selfenergies are given by |5, 

i2 1 I 111- ii'^ 11?- 11?- 



fiE,,(q )- ^ ^ y\ ^k^k-q ^ 1 ^229) 

' " 1 + |abg|\/— 2/i y ^ I iftWn — ftWk — ^k-q ihhJn + + ^k-q 2ek j ' 



k 

^2 



nEi2 (q, ^ - J2 



1 + |abg|\/— 2/i V 1^ k k k q k q yjf^^ _ fy^^ _ fiii^k-q ' +ifi'-^n + ^k + S^k-c 



and ;iE22(q, iw„) = ;iE ii(q, -za;„), ;tE i2(q, io;,,) 
/iE2i(q, ?a;„) and o^k = + (ek — A*)^ is the atomic 

spectrum in the BCS state. Note that the we have ne- 
glected the thermal Fermi factors iVk-q and A'k and that 
the corrections due to the bac kgro und interactions to the 
selfenergy (see also footnote |119| ') are approximated in 



the two-body normal state limit. This is justified in the 
BEC limit where 3> A, ksT. To the lowest order 
in perturbation theory the dressed-molecule propagator 
given by Eqs. H228|) and (|229|l can be approximated by 
the matrix 




+ eK/2 + ^ (1 - z) (i + 4^5) (1 - ^) (i + 4.;s; 



-ifkUn 



l-eK/2- 



T 



Born^ 



ihun + eK/2 + T^ 



(230) 



with 



this essentially reduces to 



T: 



Born 



(0,0,0,0) 



{4nnym^) 2a,s{B) (1 - Z)M 1 + 4 




(231) 



The derivation of Eq. H230(l follows essentially from an ex- 
pansion in |A|/|^| of the normal and anomalous molecu- 
lar selfenergies at finite energy and momentu m re tain- 
ing only the linear terms in hujn and e-K/2 jl26l |. A 
sketch of the derivation is given in Appendix A. Note 
that in Eq. (|231|l . we have made use of |Ap/4|^| = 
(iTTtf/mn,) 2acff(B)(l-Z)nn,c ^ {iTrfi^ /Tr^2a{B)n„,c, 
where aoff(i?) has been defined in Eq. 1)207(1 . Clearly, 
Eq. (|231|l implies that the molecule molecule scattering 
length in the Born approximation is given by 



2a,s{B){l-Zf h +4^1^ 



(232) 



In the BEC regime investigated experimentally until now 
for the broad resonance in ^Li, where em — —fi? /ma^{B) 



2a{B) (1 - Zy 




(233) 



However, we know that in that regime we have Z 
and ebg ^ |em|, such that 



2a{B) 



(234) 



in a very good approximation. The results of 
Eqs. (|232I233|I and H234|) has to be compared with 
the mean-field results we obtained in Eqs. H209I211() 
and (|212|1 . Moreover, it can be shown that 
the non universal correction induced by the factor 

(1 — Z)^ ^1 -I- 4^ is independent of the approxima- 
tion used in the calculation of the molecule-molecule ver- 
tex interaction. Therefore we can anticipate that the full 
dressed-molecule scattering length should obey 



am(i?) ^ 0.6 a{B) (1 - zf (^1 + j , (235) 

in the BEC limit of the crossover region after including 
higher-order corrections [sol |40| . 
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For 5 < — ebg, the approximation ei„ ~ —h^/ma?'{B) 
is not valid, and the dressed-molecules scattering length 
is expected to be 



a,„(B) ~ 0.6 acff (S) (1 ~ Zf [l 



which vanishes as len 




(236) 



at larger negative detunings 
when the wavefunction renormalization factor Z goes fast 
to 2' ~ 1. Unfortunately, the molecule-molecule scatter- 
ing length of ^Li2 molecules has not been investigated yet 
at such large negative magnetic fields [B < 650G). 

From the Bogoliubov propagator of Eq. H23()ll we find 
an energy spectrum for the dressed molecules 54] which 
is linear at low momenta 



In the remaining part of this section we want to show 
that in this regime the gas is mainly constituted by a gas 
of interacting dressed molecules. The mean- field equa- 
tion of state in Eq. (|193|l . however, needs to be modified 
after the introduction of the fluctuations that lead to the 
superfluid dressed-molecule propagator of Eq. (|230|l . In 
a self-consistent approach '14|, the total number of par- 
ticles is given by the formula 



n = 2Z?imc - 2T-r-Tr [G] + 2— -Tr [Gf] , (238) 



h(3V 



En, (K) - d{h^Ky2m^y + {h^Ky2mn,) 2TB° 



™' where the fermionic single-particle propagator Gf is given 
(237) ^jjg matrix 



\G^^^^{k,tiOn) GY^l^i^tiOn) J \ -G|,o(-k,iw„)-i y {hJ:l^ik,iiUn)+A* h^20^,iu;n) J' 

(239) 



where the fermionic selfenergies hT,lj{k,iuJn) contain the 
feedback effects of the dressed-molecules on the atoms. 
The term Tr [G] is easily evaluated according to the the- 
ory of a weak-interacting Bose gas. We have 



2Zr 



1 



2Z 



hf3V 
dq 



Tr [G] = 2Zn„ 



(240) 



(27r)3 

2Z(rw(T) + r4(r)) 



[ulicONsiEM) - vliq)NBi-EM)] 



where 



are the standard bosonic factors of the Bogoliubov trans- 
formation [s^l and n[^{T) is the noncondensed density 
of dressed molecules at temperature T. The density of 
Eq. H240|l corresponds to the bare-molecular contribution 
to the total density. It is very small because in the BEC 
limit that we are considering in a broad resonance, we 
have Z ~ 0. therefore, we show that the main contri- 
bution in the equation of state in Eq. 1)2381) comes from 
Tr [Gf ] . To calculate this trace we have first to calculated 
the fermionic selfenergies. Upon neglecting higher order 
contributions, we ultimately find |33| 



<(q) = <(q) - 1 



q 

2m„ 



q^Boin 



in "Tiic 

2-Bm(q) 



^^m(q) 



(241) 




where 



Note that we have neglected the off-diagonal fermionic 
selfenergies because in the limit under consideration one 
(243) can show that ^^^2(0,0) ~ |A|(|A|V2|Aip)- From 
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Eq. I|242|) we get the expression for Gf^i in the BEC 
hmit ,33i| 



Gtai(k,LJs) 



ihujs + /i — Ek — 



|AP+A'^ 



(244) 



where we have neglected a term of order A' with 
respect to Note that Eq. 1244|) has the same formal 
structure of the corresponding BCS expression with 
the replacement (A^ + A'^). Accordingly, the trace 

of the fermionic propagator can be rewritten in the BCS- 
like form 



1 



h(3V 



2 Tr [Gf] 



(245) 



(27r) 



d^k 1 



(ek - f^y 



Expanding the right-hand side of this equation as we did 
for the mean-field BCS equation of state in Eq. H196|) we 
have 



1 



h(3V 



2 Tr [Gf] ~ (|A|2 + A'- 



(246) 



Moreover, using 2/i ~ |e„i| and the definition of two-body 
Z we find 



1 



TV[Gf] 



2Z (rime + »^m') V 

2 (1 - Z) (n^e + n'J . 



(247) 



Joining the two contributions of Eqs. H247|l and H24()|l we 
have 



n~2(n„,c(r)-f <(T)). 



(248) 



This result holds asymptotically for ksT <C |em| and it 
represents the extension to a two-channel model in the 
case of a broad resonance, of the description of the BEC 
limit in Ref. jsjjs^ based on the composite bosons of the 
single- channel model. However, it is important to stress 
that these results are only valid in the lowest order of 
perturbation theory. Furthermore, we use a definition of 
the probability Z that is only true in the asymptotic BEC 
limit. The calculation of Z in the strong-coupling region 
requires a more general approach as shown in Ref. 54] . 

Nevertheless, our analysis of the BEC limit, shows the 
way to connect the two models in the limit of a very broad 



resonance, where the single-channel model accounts fairly 
well for the thermodynamics of the gas inside the BEC- 
BCS crossover region. 



VI. CONCLUSIONS AND OUTLOOK. 

In this paper we have shown that the dressed-molecule 
method represents an effective and consistent approach 
in order to describe the physics of ultracold Fermi gases 
near a Feshbach resonance. The dressed molecule is 
the real physical entity of the BEC-BCS crossover in 
atomic Fermi gases near a Feshbach resonance, because 
the wavefunction of the Feshbach molecules is always a 
coherent superposition. In this formulation the informa- 
tion about the mixing of the two-channels is preserved 
at each stage in the many-body calculation. Therefore it 
has the advantage that it can be applied also to medium 
and narrow reconances, when i]^ < ep, where the single- 
channel approximation is expected to fail. However, a 
proper treatment of the narrow resonance case requires 
the inclusion of the finite range corrections in the theory 
This will be addressed elsewhere. 
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Appendix A. 

In this appendix we present the detailed calcula- 
tion which leads to the dressed-molecule propagator of 
Eq. (|230|l from Eq. (|228|l . In order to make the derivation 
more transparent we discuss first the case when abg = 0, 
because the inclusion of the background interactions does 
not change qualitatively the discussion. Neglecting the 
background scattering length corrections, the molecular 
selfenergies of Eq. (|229|l reduce to 
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«s„«,,.:h.)4 T"-", — +^ . (249) 



k 

^2 



Vv- / 1 , 1 



?li;i2(q,lW„) = — ^ <Uk?^kUk-qWk-q ( — 



I 

where, in the definition of the gap of Eq. ()188|l entering with the coefRcients of the expansion given by 

the definition of the coherence factors Vk and uu, the 

contribution due to the background interaction has to be 

neglected as well. In the BEC limit, where 2/i ~ em ^ 

ksT, A, we can expand the selfenergies of Eq. 1249|) for 

small values of |q| and uJn as 

fiEii(q,iw„) =Ao + Ai ihuJn + A2 eq/2 + 0(w2jq|4) 
;iSi2(q,ic^„) =Bo + OiLul\q\^). (250) 
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° y 2^2ek 2c.k 1^ ^4c.3' y 4 ^3 

I[ (^u-M)(2(ek-Mf-|AP) 2.k|AP(8(ek-Mf + 3|AP) | ^2 |^| 

^^"^n 44 ^ ^^^^i /' ""'-y^^^. ^'''^ 

Substituting these results in Eq. (|228|) and using the gap equation of Eq. IjigOf) the dressed-molecule propagator can 
be rewritten as 



-ihLo,, + eK/2 - ^ Ek + + M eK/2 ^ Ek ^ 



f Ek IS ^^'^^ + eK/2 - % Ek IS - ^1*^^" + A2 eK/2 

k k 

(252) 

This expression can be simplified further by noting that in the BEC limit one has 

V ^i^i \ Z ) 4|mI ' ' \ Z )• ' ' 

in the leading order of |Ap/|/i|. Therefore, we get 



, 1 /-iftw„ + eK/2+#T(l-^) ^-Pri^-Z) 
^(1-^) + eK/2 + ^ 

r 



^(1-^) zfi^„ + eK/2 + ^(l-Z), 



This derivation can be generalized rather straight- 
forwardly in order to include the effects of the back- 
ground scattering length in the molecular selfenergies of 
Eq. if^ . The coefficients of Eq. differ only by 

an overal factor 1/ (l -I- |abg|\/— 2/7) , and Z, em, A and 
the gap equation, must be replaced everywhere with the 
more general expressions that include the background 



scattering length corrections. Note, however, that the 
background scattering corrections in the selfenergies of 
Eq. (|229|l are valid only in the BEC hmit. Therefore, in 
order to semplify the dressed-molecule propagator as we 
have done in Eq. (|252(l . we cannot use the general gap 
equation of Eq. p9UI) . Alternatively, we expand the co- 
efficients of Eq. H251|) in powers of A/|/i| and then we use 
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the expanded version of the gap equation of Eq. H195|) . 

Ultimately we find 

^ 1 /-ifii^„ + eK/2 + ^(l-Z)(l + 4^ 



1 + 4, 



^(l-Z)(l + 4v^) 



(255) 



as given 



in Eq. if^ . 
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